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BRIEF  SUMMARY 
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The  purpose  of  this  experiment  was  to  verify  the  analysis  of  radial 
equilibrium  for  compressible  flow  through  a stationary  cascade,  where 
radial  equilibrium  is  defined  as  uniform  axial  mass  flow.  Experimental 
determination  of  the  effect  of  compressibility  on  analytical  predictions 
of  flow  angle  was  of  primary  importance.  Previous  analytical  work  is 
sumnarized  and  those  results  shown. 

A specially  designed  annular  flow  test  apparatus  was  constructed  for 
this  work  and  its  features  are  described  in  detail. 

Finally,  the  results  of  detailed  flow  field  surveys  using  three 
different  blade  angle  settings  are  summarized  and  compared  with  both  the 
free  vortex  and  the  Engelman  analysis.  In  order  to  aid  visualization  of 
the  real  flow  field  as  measured,  various  computerized  graphical  presentations 
of  the  data  are  shown.  In  addition,  these  survey  results  are  compared  with 
analytical  predictions  from  a two-dimensional  computer  program  developed  by 
the  General  Electric  Company  and  a three-dimensional  computer  program  developed 
by  AiResearch  Manufacturing  Company  of  Arizona. 
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CHAPTER  I 


INTRODUCTION 

The  high  output  of  turbine  engines  depends  on  many  factors.  One  of 
the  more  significant  ones  is  predicting  the  flow  characteristics  through 
an  annular  cascade.  Any  improvement  in  accuracy  of  the  methods  available 
would  contribute  greatly  to  the  design  of  such  engines. 

Traditionally,  in  the  early  1900's,  the  design  of  compressor  or 
turbine  stages  was  based  on  the  pitch-line  calculation  using  the  mean 
camberline  as  the  reference  for  the  inlet  and  exit  angles.  Later,  as 
the  rotational  speed  and  loading  per  blade  increased,  two  additional 
important  factors  were  included  in  the  design.  These  are  the  influence 
of  the  centripetal  force  and  the  effect  of  the  finite  thickness  of  the 
blade  airfoil  . 

The  effect  of  high  rotational  speed  is  to  cause  the  flow  to  gravi- 
tate toward  the  shroud  area  of  the  blade,  thus,  starving  the  hub  area  of 
the  blade.  The  net  result  is  the  undesirable  effect  of  not  fully  utiliz- 
ing the  entire  span  of  the  blade.  To  overcome  this  maldistribution,  the 
concept  of  designing  the  blade  profile  to  reach  a radial  equilibrium  was 
evolved  in  the  1 930 ’ s . Radial  equil ibrium  means  that  no  radial  flow 
exists  in  the  flowpath  and,  therefore,  there  exists  uniform  axial  mass 
flow.  The  basic  idea  is  simply  to  "underexpand"  the  flow  near  the  tip 
section  and  to  "overexpand"  near  the  hub  section.  This  controlled 
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expansion  creates  an  increasing  pressure  distribution,  along  the  radial 
direction,  which  is  utilized  to  maintain  uniformity  of  the  mass  flow  in 
the  axial  direction.  In  this  manner,  the  loading  of  the  blade  in  the 
span  direction  is  kept  uniform,  at  least  from  a design  viewpoint.  The 
calculation  of  the  necessary  twist  of  the  blade  along  the  spanline  direc- 
tion is  designated,  generally  speaking,  the  radial  equilibrium  theory. 

(11* 

Among  those  who  have  contributed  greatly  are  Wu  and  Wolfensteinv  ' , 

Wu^,  Wu^,  Stanitz^,  Wu  and  Brown^,  Huppert  and  MacGregor^, 

Smith  and  Traugott^,  Shepherd^,  Giamati  and  Finger^,  Mereiros  and 

Hood^),  Herzig  and  Hansen^^,  Smith^^,  Hawthorne  and  Horlock^^, 
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and  Hawthorne  and  Ringrose'  . These  are  but  a few  of  the  possible 
citations,  and  Marble^15,  pp'  develops  the  entire  three-dimen- 

sional flow  treatment  in  great  detail. 

A second  advancement  in  the  design  procedure  is  the  use  of  a mapping 
technique  in  conjunction  with  the  potential  flow  theory  to  calculate  the 
streamlines  through  a cascade.  The  turning  of  flow  can  thus  be  more 
realistically  assessed  with  the  airfoil  thickness  taken  into  account.  A 
drawback  is  that  such  a flow  calculation  is  restricted  to  a two-dimen- 
sioral  or  linear  cascade.  For  an  annular  cascade,  however,  there  has 
been,  up  to  the  present,  no  theory  advanced  to  account  for  the  curvature 
effect.  Consequently,  a compromise  procedure  is  to  calculate  the  flow 
characteristics  at  various  radial  lines  for  an  annular  cascade.  The 
resultant  flow  distributions  are  stacked  upon  each  other  to  form,  at 
least  as  a first  order  solution,  a composite  picture  for  an  annular 
cascade. 

♦Numbers  in  parentheses  ( N , p . M ) refers  to  item  number  (N)  in  the  bibli- 
ography and  the  specific  page  (M)  of  that  reference  as  appropriate 


More  recently,  attempts  have  been  made  by  Dodge^^  to  use  a numer- 


al ical  scheme  to  calculate  the  flow  distribution  in  a channel  formed  by 

two  adjacent  blades  and  an  upper  wall  representing  the  shroud  and  a 
lower  wall  representing  the  hub.  This  numerical  scheme  takes  into 
account  the  effect  of  viscosity,  wall  roughness,  etc.,  and  its  development 
is  still  in  progress. 

The  performance  improvement  due  to  the  present  design  capability  to 
include  all  these  considerations  for  both  compressor  and  turbines  stages 
has  been  demonstrated,  in  a general  sense,  by  increased  propulsive  effi- 
ciency. For  gas  turbine  units  of  small  size  such  as  Auxiliary  Power  Units 
(APU),  these  modern  tools  are  yet  to  be  fully  incorporated  into  their 
V design.  This  is  due  partly  to  a size  limitation  and  partly  to  the  greatly 

increased  cost  versus  possible  gains  in  their  efficiency.  For  the  reasons 
of  ascertaining  the  usefulness  of  the  concepts  of  radial  equilibrium  the 
two-dimensional  cascade  flow  theory,  and  tr.e  three-dimensional  viscous 
calculations,  the  present  investigation  was  undertaken  in  part  to  deter- 
mine the  extent  to  which  they  can  or  should  be  applied  to  the  design  of 
small  units.  For  economical  reasons,  the  blades  of  the  annular  cascade 
in  this  investigation  were  chosen  straight  so  as  to  ascertain  the  flow 
departure  from  that  of  radial  equilibrium.  In  addition,  a major  goal  of 
this  investigation  was  to  verify  Engelman's  analysis  of  radial  equilibrium 
in  compressible  flow. 

It  has  been  concluded  by  several  authors.  Binder^’  p'  and 

Shepherd^8,  p>  205,  that  radial  equilibrium  is  a desirable  condition 

to  achieve  with  the  machine  and  blade  designs. 

Using  the  proper  assumptions  to  simplify  the  basic  mathematical 
description  of  the  physical  situation,  the  calculation  of  radial 
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distribution  of  flow  can  be  reduced  to  one  basic  equation: 


As  is  well  known,  e.g.  see  Velkoff^^’  , the  simplest  and 

therefore  most  used  solution  scheme  is  to  make  the  assumption  that 
between  blade  rows,  the  actual  flow  can  be  represented  by  a free  vortex 
flow.  However,  such  an  assumed  flow  field  will  not  always  produce  radial 
equilibrium  and  Professor  H.  W.  Engelman'  ' has  suggested  that  two  solu- 
tions are  appropriate;  i.e.,  one  for  incompressible  flow  and  another  for 
compressible  flow.  In  his  derivation  and  proposed  solutions,  the  achieve- 
ment of  radial  equilibrium  is  the  primary  condition  and  the  assumptions 
made  are: 

(1)  constant  enthalpy 

(2)  constant  entropy 

(3)  mass  is  conserved 

The  only  difference  between  the  two  solutions  is  that  in  one  case 
the  fluid  is  assumed  to  be  incompressible,  while  in  the  other  the  fluid 
is  considered  as  an  ideal  compressible  gas. 

A dimensionless  solution  was  derived  by  Engelman  and  numerical  solu- 
tion of  the  resulting  equations  was  carried  out  by  Bogus  . The  result 
for  the  incompressible  fluid  is  a free  vortex,  while  the  result  for  the 
compressible  fluid  is  very  different  in  the  transonic  and  supersonic 
regimes.  A final  observation  is  that  the  compressible  result  is  very 
close  to  a free  vortex  in  the  low  subsonic  range. 

Complete  derivations  of  the  basic  equation,  and  the  numerical  solu- 
tions are  included  in  Appendix  A. 
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In  an  attempt  to  obtain  experimental  verification  of  these  analytical 
results,  a small  rig  was  fabricated  and  tested  at  Robinson  Laboratory 
by  Bogus  in  1970.  Due  to  the  channel  configuration  and  size  relative 
to  the  available  measurement  system,  the  results  were  inconclusive. 

In  summary,  the  primary  purpose  of  this  investigation  was  to  obtain 
sufficient  experimental  data  to  either  prove  or  disprove  the  validity  of 
the  Engelman  radial  equilibrium  analysis. 
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CHAPTER  II 


EXPERIMENT 


Design  of  Test  Apparatus 

Since  no  existing  test  configuration  was  available  which  could  be 
used  for  this  purpose,  an  entirely  new  apparatus  had  to  be  designed  and 
built.  As  a starting  point,  it  was  clear  from  the  previous  attempt  by 
Bogus^)  that  at  least  a 12-inch  outer  diameter  annulus  would  be 
required.  This  was  considered  necessary  to  prevent  probe  and  endwall 
effects  from  dominating  the  results.  Such  a large  channel  would,  of 
course,  require  a significant  increase  in  facility  pumping  capacity  com- 
pared to  the  Ohio  State  University  capability.  The  test  facilities  of 
the  Air  Force  Aero  Propulsion  Laboratory  at  Wright-Patterson  Air  Force 
Base  include  altitude  chambers  capable  of  handling  the  through  flow  of 
J-85  class  engines  up  to  approximately  50,000  feet.  Thus  the  operation 
of  an  annular  flow  cascade  in  this  class  was  considered  feasible  and  a 
preliminary  design  layout  was  made. 

The  parameters  of  interest  which  were  to  be  measured  at  various 
points  in  the  channel  were  total  and  static  pressures,  total  temperature, 
and  flow  angle.  To  measure  these,  a wedge  probe  was  chosen  and  was  used 
for  all  measurements.  This  probe  provides  three  pressure  taps;  one  total 
pressure  and  two  static  pressures,  one  in  each  face  of  the  wedge.  The 
flow  angle  was  measured  by  rotating  the  probe  until  the  two  wedge  face 


tap  pressures  balanced  each  other.  A checkout  setup  was  fabricated  to 
assure  the  accuracy  of  this  probe  for  use  in  measuring  flow  angle.  The 
details  of  this  checkout  are  described  in  Appendix  B. 

An  atmospheric  inlet  was  chosen;  its  design  therefore  required  some 
form  of  bellmouth  to  assure  a uniform  axial  flow  field.  Fortuitously, 
there  existed  a surplus  J-85  engine  test  bellmouth  with  a 16-inch  inside 
diameter  and  preliminary  calculations  indicated  that,  with  a reasonable 
hub-tip  ratio,  the  laboratory  air  supply  system  could  handle  the  result- 
ing flow.  Therefore,  the  preliminary  layout  was  based  on  a 16-inch 
channel  outer  diameter  with  an  11 -inch  centerbody  outer  diameter  thus 
giving  a 2 1/2-inch  blade  height  and  a 0.6875  hub  to  tip  ratio. 

In  order  to  produce  the  whirling  flow  field  required  for  this  exper- 
iment, nontapered,  nontwisted  blades  were  desired.  This  would  produce 
a flow  field  which  was  not  specifically  tailored  to  fit  a certain  radial 
equilibrium  model  or  theory,  but  rather  one  which  could  then  be  surveyed 
to  ascertain  its  characteristics.  In  addition,  of  course,  the  Engelman 
analysis  predicts  that  such  blading  will  produce  uniform  axial  mass  flow 
distribution  at  high  Mach  numbers  and  verification  would  then  be  a 
straightforward  experiment.  After  giving  brief  consideration  to  design- 
ing and  fabricating  a special  blade,  it  was  concluded  that  surely  some 
existing  blade  could  be  found  which  would  perform  adequately.  As  a 
result  of  work  sponsored  at  Teledyne  CAE^  ’ ^ by  the  U.  S.  Army,  an 
extruded  aluminum  blade  design  existed  which  was  considered  quite  satis- 
factory for  the  intended  purpose.  A closeup  of  this  contour  is  shown  in 
Figure  1 and  the  profile  coordinates  are  given  in  Figure  2.  Permission 
was  obtained  from  the  Army  to  use  the  extrusion  die  to  produce  enough 
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Figure  2 Vane  Coordinates 
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strip  stock  for  three  rows  of  blades.  The  design  point  blade  angle  for 
this  profile  was  69.7°  and  one  blade  row  was  set  at  this  angle.  Two 
more  rows  were  fabricated  with  settings  of  65°  and  60°  to  provide  options 
for  a broad  range  of  test  conditions.  Pictures  of  the  blade  rows  are 
shown  as  Figures  3,  4,  5,  and  6 with  engineering  drawings  shown  in 
Figures  7 through  12.  Tables  1,  2,  and  3 show  the  geometric  and  aerody- 
namic design  data. 

Having  established  the  channel  cross  sectional  configuration,  only 
channel  length  remained  to  be  set.  The  goal  was  at  least  24  inches  with 
as  many  axial  measurement  stations  as  practical.  As  shown  in  Figure  13, 
the  final  as-built  channel  length  was  approximately  twenty-eight  inches 
and  eleven  measurement  stations  were  incorporated  at  two  inch  centerline- 
to-centerl ine  spacing.  Figure  14  shows  the  numbering  system  used  for 
stations  and  centerbody  static  pressure  taps. 

Also  shown  in  Figure  13  are  the  remaining  construction  details  of 
the  assembly,  including  two  unique  and  patented  features.  The  first  was 
devised  to  permit  circumferential  traversing  of  the  probe  with  minimum 
aerodynamic  disturbance  and  the  second  relates  to  the  automatic  wrenchless 
coupling/decoupling  of  centerbody  static  pressure  tap  readout  lines. 

In  order  to  move  the  probe  circumferentially,  slots  would,  of  course, 
be  required  in  the  channel  wall.  Cutting  such  slots  would  seriously 
interrupt  the  channel  boundary  and  create  undesirable  leakage  flows.  So 
after  making  many  preliminary  sketches  and  reviewing  conceptual  ideas,  it 
was  decided  to  make  a liner  for  the  rig  which  would  serve  as  the  channel 
outer  wall  and  which  could  rotate  inside  the  slotted  outer  structure. 

This  can  be  seen  in  Figures  15  and  16.  In  this  way,  only  simple  radial 
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Fiqure  7 69.7°  Blade  Row  Detail 


Parameter 

Percent  Span  at 

Exit 

0° 

50" 

100" 

Inlet  Angle  (degrees) 

0.0 

0.0 

0.0 

Exit  Angle  (degrees) 
(Pressure  Side) 

69. 7‘ 

69.7° 

69.7C 

Axial  Chord  (in) 

.830 

.830 

.830 

Spacing  (in) 

0.72 

0.88 

1 .05 

Throat  (in) 

.2126 

.4625 

.5875 

Trailing  Edge  Radius  (in) 

0.015 

0.015 

0.015 

Leading  Edge  Radius  (in) 

0.070 

0.070 

0.070 

Aspect  Ratio  (Span/Chord) 

-- 

0.576 

Spacing/Chord  Ratio 

0.87 

1 .06 

1.27 

Number  of  Vanes 

-- 

48 

-- 

k 


Table  2 


Summary  of  Vane  Aerodynamic 
and  Geometric  Design  Data 

65°  Blades 


Parameter 

3ercent  Span  at 

Exit 

0% 

50? 

S'* 

O 

O 

Inlet  Angle  (degrees) 

4.9° 

4.9° 

4.9° 

Exit  Angle  (degrees) 

65° 

65° 

65° 

Axial  Chord  (in) 

.905 

.905 

.905 

Spacing  (in) 

0.691 

.85 

1 .01 

Throat  (in) 

0.266 

0.356 

0.438 

Trailing  Edge  Radius  (in) 

0.015 

0.015 

0.015 

Leading  Edge  Radius  (in) 

0.070 

0.070 

0.070 

Aspect  Ratio 

0.576 

-- 

Spacing/Chord  Ratio 

0.764 

0.939 

1.116 

Number  of  Vanes 

■■  ■ 

50 

4 

CO  WTO  UR  AT  11,00  O' A. 


Figure  12  60°  Blade  Row  Sections 


Table  3 


Summary  of  Vane  Aerodynamic 
and  Geometric  Design  Data 

60°  Blades 


l 


Percent  Span  at  Exit 


Inlet  Angle  (degrees) 


Exit  Angle  (dearees) 
(Pressure  Side) 

Axial  Chord  (in) 


Spacing  (in) 


Throat  (in) 


Trailing  Edge  Radius  (in) 


Leading  Edge  Radius  (in) 


Aspect  Ratio  (Span  Chord) 


Spacing/Chord  Ratio 


Number  of  Vanes 
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Figure  15  Outer  Shell  With  Wood  Liner  Removed 


holes  were  required  in  the  liner  to  permit  probe  penetration  and  those 
not  in  use  could  be  plugged  with  simple,  nonflow-interfering  slugs.  This 
design  concept  was  awarded  United  States  Patent  No.  4,031,749  on 
28  June  1977.  The  liner.  Figure  16,  was  made  of  laminated  mahogany  for 
dimensional  stability  and  was  approximately  three  and  one  half  inches 
thick  to  permit  flexibility  in  channel  wall  contour,  if  later  required. 
The  penetrations  through  the  liner  consisted  of  metal  inserts  in  the  wood 
which  were  drilled  large  enough  to  provide  clearance  for  the  probes. 

In  addition,  the  flow  channel  side  was  tapped  for  either  a probe  support/ 
seal  bushing  or  a solid  cylindrical  plug.  These  were  made  of  teflon  to 
reduce  friction  on  the  probe,  and  a typical  one  is  shown  in  Figure  17. 

The  liner  was  sealed  to  the  outer  structure  at  each  end  with  an  0-ring 
made  of  hollow  plastic  tubing. 

With  regard  to  the  wrenchless  static  pressure  line  connections, 
this  was  achieved  by  putting  an  elastomeric  0-ring  around  each  line  to 
act  as  a static  seal  as  the  line  crossed  a bulkhead  in  the  rear  center- 
body  support  assembly.  This  bulkhead  is  shown  in  Figure  18.  To  allow 
for  any  slight  misalignment  due  to  fabrication  tolerance  or  temnerature 
changes  during  a test  run,  each  0-ring  was  contained  in  a spring-loaded 
carrier  which  permitted  slight  axial  and  radial  displacements.  This 
connection  concept  was  awarded  United  States  Patent  No.  4,111,463  on 
5 September  1978. 

The  entire  assembly  was  designed  for  easy  teardown  and  quick  access 
to  the  flow  channel  and  all  instrumentation.  Standard  tooling  toggle 
clamps  were  used  to  retain  both  the  bellmouth  and  the  circumferential 
saddle.  Thus,  without  any  wrenches,  the  axial  saddle  station  could  be 
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changed  and  the  inlet  bellmouth  assembly  removed.  Once  the  bellmouth 
assembly  was  off  the  main  assembly,  the  blade  row  could  be  slipped  out, 
again  without  tools.  This  gave  access  to  the  channel  and  probe.  Then, 
with  only  the  removal  of  one  hex  nut,  the  entire  centerbody  assembly 
could  be  removed.  As  described  above,  the  static  pressure  lines  were 
automatically  broken  at  the  rear  of  the  centerbody  without  any  require- 
ment for  tools.  These  details  are  shown  in  Figures  19  through  26. 

The  centerbody  static  pressure  line  installation  is  shown  in 

Figure  24.  All  wall  static  pressure  taps  used  in  this  experiment  were 
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machined  in  accordance  with  the  recommendations  of  Rayle'  as  reported 
by  Dean^4,  pp‘  ® . The  centerbody  taps  are  shown  in  Figure  27. 

The  first  thirteen  are  spaced  one-half  inch  apart  and  the  others  are 
spaced  one  inch  apart.  A closeup  of  a typical  outer  wall  static  tap  is 
shown  in  Figure  28.  The  teflon  bushing  adapts  the  metal  pressure  tap  to 
the  liner  penetration  ports  in  a manner  similar  to  that  used  for  the 
wedge  probe,  and  these  are  shown  as  installed  in  Figure  29. 

The  circumferential  saddle  assembly  shown  in  Figures  30  and  31  was 
designed  to  permit  mounting  two  probe  actuators  manufactured  by  Northern 
Research  Company.  These  actuators  are  hydraulic  servo  positioning 
mechanisms  which  can  be  remotely  set  from  the  test  cell  control  room. 
They  accurately  and  repeatably  set  the  radial  and  angular  positions  of 
any  probe  inserted  in  their  collet  assembly.  Details  of  the  calibration 
procedures  used  on  the  actuator  are  contained  in  Appendix  C.  A three- 
channel  system  was  purchased  since  no  similar  equipment  existed  in  the 
Laboratory,  and  the  flexibility  for  future  use  was  considered  worth  the 
extra  cost.  In  the  final  arrangement  for  this  experiment,  only  one 
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actuator  was  used.  The  system  consists  of  the  actuators,  a hydraulic 
power  supply,  and  a remote  control /readout  panel.  These  are  shown  in  a 
bench  functional  checkout  setup  in  Figure  32.  The  actuator  is  shown  as 
mounted  on  the  circumferential  saddle  assembly  on  the  rig  in  Figures  31 
and  33.  A complete  numerical  index  of  all  rig  fabrication  drawings  is 
shown  in  Appendix  D. 

Experimental  Setup 

The  entire  test  assembly  with  probe  installed  is  shown  in  Figures  34 
and  35  as  set  up  in  Cell  24  of  Building  18E  at  W^ight-Patterson  Air  Force 
Base,  Area  B.  A conical  transition  piece  was  fabricated  to  connect  the 
test  apparatus  to  an  existing  altitude  engine  stand  via  approximately  18 
feet  of  16-inch  Schedule  40  pipe.  All  rig  controls  and  instrumentation 
readout  devices  were  remote  mounted  in  the  cell  control  room.  The  only 
parameter  which  had  to  be  read  directly  from  the  rig  itself  was  the  cir- 
cumferential position  of  the  saddle.  This  scale  was  graduated  in  degrees 
with  zero  being  the  vertical  centerline  of  the  saddle.  Therefore,  since 
the  saddle  contains  two  actuator  pads,  the  actual  circumferential  position 
of  each  probe  radial  line  of  action  was  offset  from  zero.  These  pads 
were  machined  as  close  together  as  possible;  but  due  to  the  actuator  size, 
each  had  to  be  offset  7 1/2  degrees;  i.e.,  there  was  a 15-degree  included 
angle  between  probe  radial  paths.  These  details  are  shown  in  Figure  13. 
The  relative  location  of  the  probe  and  blade  row  with  probe  installed 
at  station  one  is  described  in  Figures  36  and  37. 

Except  for  total  temperature  readout  from  the  probe,  the  control 
room  arrangement  is  shown  in  Figure  38,  while  Figure  39  shows  the  Brown 
recorder  used  for  temperature  readout  and  recording.  All  probe  pressures 
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Figure  35  Test  Assembly  Installation  - Left  Side 


Figure  36  PROBE  - CHANNEL 
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Figure  37  Spatial  Relation  of  Stations  to  Bladinq,  and  Geometric  Parameters  for  the  3 Blade  Settings  Used 


were  read  on  Kollsman  aneroid  gages  graduated  in  inches  of  water  which 
were  calibrated  by  the  Base  PME  (Precision  Measuring  Equipment)  Shop. 

The  remaining  gages  on  the  panel  were  used  for  outer  channel  wall  static 
readings,  while  the  inner  wall  static  taps  were  connected  to  a manometer 
board  also  shown  in  Figure  38. 

The  probe  angular  and  radial  (linear)  positions  were  read  from  a 
digital  voltmeter  which  was  part  of  the  actuator  package.  Using  the 
calibration  table  in  Appendix  C,  the  actual  geometrical  position  could  be 
determined  or  set. 

Da ta  Acquisition  Technique 

The  raw  data  taken  at  each  circumferential  position  consisted  of: 
Wedge  Probe 

Ptotal 

/ P = P - P 

total  static 

"'"total 
Flow  Angle 
Radial  Position 

Wall  Static  Pressures  (both  inner  and  outer) 

Chamber  (Receiver)  Pressure 
Cell  (Air  Inlet)  Temperature 

As  stated  before,  it  was  not  considered  cost  or  time  effective  to 
install  and  checkout  the  required  data  automation  equipment  on  such  a 
relatively  simple,  steady-state  experiment.  Therefore,  all  data  were 
recorded  by  hand  on  prepared  data  sheets,  and  two  sample  complete  runs 
are  shown  in  Appendix  E.  These  sample  data  sheets  show  the  total  scheme 
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used  for  recording  all  data  taken.  Note  that,  on  the  data  sheets  P1 ; 
P^  and  P^  are  the  wedge  face  taps. 

The  standard  nomenclature  for  probe  circumferential  location  used 
for  all  data  taken  during  this  experiment  is  shown  in  Figure  40  and  is 
described  below. 

As  shown  in  Figures  7,  9,  and  11,  each  of  the  three  cascades  was 
located  by  a dowel  pin  such  that  a trailing  edge  was  always  on  the  verti- 
cal centerline  of  the  test  assembly.  In  addition,  the  circumferential 
"saddle"  zero  reference  was  located  on  the  vertical  centerline.  Now, 
since  the  saddle  was  designed  and  machined  to  permit  simultaneous  mount- 
ing of  two  probe  actuators,  the  included  angle  between  probe  lines  of 
action  was  15  degrees.  Therefore,  when  the  saddle  circumferential  posi- 
tion indicating  pointer  is  set  at  zero,  the  actual  Drobe  positions  are 
7 1/2°  right  or  left  from  zero.  Finally,  right  and  left  were  established 
as  seen  when  facing  the  cascade  inlet.  So,  during  data  taking,  the  cir- 
cumferential position  was  always  read  directly  from  the  scale  and  a 
proper  allowance  was  made  for  actual  probe  locations  during  data  reduc- 
tion. For  all  runs  conducted,  the  wedge  probe  actuator  was  mounted  on 
the  right  saddle  pad  so  that  the  following  conversion  was  apolicable  to 
the  raw  data  sheets:  0°  indicated  7 1/2°  right. 

The  procedure  used  during  a run  at  a given  circumferential  position 
was  to  set  the  desired  radial  traverse  position  and  then  rotate  the 
wedge  probe  until  the  pressure  difference  between  the  two  side  static 
taps  was  zero.  At  this  point,  the  flow  angle  was  recorded  and  all  other 
readings  taken  in  sequence.  The  standard  radial  traverse  positions  used 
are  shown  in  Table  4.  A "short"  traverse  was  also  developed  later  as 
data  became  available,  and  comparative  plots  made  which  showed  little 
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Figure  40  Probe  Circumferential  Location  Nomenclature 
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Table  4 


— 

Linear 

Vol tage 

Position 

Reading 

0.29 

1 .956 

0.39 

2.193 

0.49 

2.444 

0.59 

2.701 

0.69 

2.966 

0.79 

3.218 

0.89 

3.448 

0.99 

3.686 

1 .09 

3.916 

1.19 

4.160 

1.29 

4.39S 

1 .39 

4.642 

1 .49 

4.894 

1.59 

5.133 

1 .69 

5.378 

1.79 

5.625 

1 .89 

5.866 

1 .99 

6.113 

2.09 

6.358 

2.19 

6.627 

2.29 

6.857 

2.39 

7.066 

2.49 

7.227 
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need  for  the  close  spacing  in  the  center  of  the  flow  channel.  A sample 
of  these  plots  is  shown  in  Appendix  F. 

During  the  trial  runs,  it  was  observed  that  total  pressure  readings 
across  the  channel  were  not  uniform  and  that  recirculation  was  apparently 
occurring  inside  the  channel  due  to  exit  conditions.  Considerable  effort 
was  then  expended  to  eliminate  this  characteristic  with  the  final  result 
being  a "backflow  plate"  or  flow  fence  installed  at  the  outer  diameter 
of  the  channel  exit.  The  proper  plate  bore  for  each  blade  row  was  deter- 
mined by  an  experimental  trial -and-error  procedure.  Plots  of  total 
pressure  and  flow  angle  vs  probe  position  as  a function  of  backflow 
plate  bore  are  shown  in  Appendix  G.  As  can  be  seen  from  these  data, 
a slight  fence  height  change  makes  a very  significant  difference  in  the 
pressure  profile.  For  the  60°  blade  cascade,  no  plate  was  required. 

Once  this  was  settled,  actual  traverse  data  were  taken  over  a wide 
range  of  receiver  pressures  (flow  rates)  and  at  many  axial  and  circumfer- 
ential positions  in  the  channel.  All  three  blade  cascades  were  even- 
tually tested  and  the  results  are  presented  in  Chapter  III. 


59 


CHAPTER  III 


RESULTS 

The  experimental  results  obtained  are  summarized  and  discussed  in 
this  section.  An  analysis  of  experimental  error  associated  with  the 
instrumentation  and  procedures  used  is  presented  in  Appendix  L.  Data  was 
taken  for  all  three  blade  angle  settings:  60°,  65°,  and  69.7°.  Computer- 
generated plots  of  the  basic  data  are  shown  for  selected  representative 
runs  in  Figures  41-43,  45-47,  49-51,  53-55,  57-59,  and  61-63.  The  curve 
fitting  method  was  a spline  technique.  To  serve  as  a guide  to  the  range 
of  profiles  covered,  the  following  tabulation  organizes  the  entire  set 
of  six  samples  with  respect  to  probe  circumferential  location.  Recall 
that  the  circumferential  location  is  specified  in  degrees  right  or  left 
when  facing  the  inlet  to  the  cascade  and  that  the  zero  reference  is  the 
cascade  vertical  centerline  which  always  has  a blade  trailing  edge 
indexed  on  it  by  means  of  a dowel  pin. 


Table  5 

Circumferential  Location  of  Sample  Traverses 


Blade 

Station  No 

Angle 
(degrees ) 

1 

3 

69.7 

2 . 5°L 

0° 

65 

1 0°R 

7 . 5°R 

60 

2 . 5°R 

5 . 5°R 
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PROBE  POSITION  UN  PROM  CENTERBODY) 
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Blade  Rngle  69.7° 
Station  1 
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1 

Blade  Angle  69.7°  j 

Station  1 


Figure  42  aP  VS  Probe  Position 
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PROBE  POSITION  (IN  FROM  CENTERBODT ) 


PLOW/BREfl  (LB/SEC/FT**2) 

Figure  44  Flow  VS  Probe  Position 
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Blade  Angle  65.0° 
Station  1 
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Blade  Angle  65.0° 
Station  1 


Blade  Angle  65.0° 
Station  1 
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Figure  47  Flow  Angle  VS  Probe  Position 


68 


r~> 

X X) 

X 

VP 

UJ 

IP 

X 

<X3 

X 

X 

X 

UJ 

X 

© 

if. 

ft 

p 

X 

CD 

o 

P 

X 

*-* 

X fs. 

f"  > 

O' 

«-« 

~ > 

X 

X 

X 

UJ 

a 

X 

X 

X 

X 

V' 

X 

X 

CO 

X 

*V 

IP 

* ft 

<o 

X 

cr 

O' 

J 

IP 

cr 

X 

X 

X 

X 

X 

-3 

p 

IP 

© 

O' 

ft 

X. 

ft 

x p 

vD 

a. 

CD 

-j 

UJ  P 

o 

X 

X 

<D 

OJ 

a 

X 

X 

J 

ft 

X 

p 

X 

CD 

in 

cr 

rr  cr 

cr 

O' 

CD 

o 

CD 

r— 

D 

o 

CD 

CD 

o 

X 

X 

X 

ft 

■T 

n 

ft) 

o 

X 

a 

p 

• 

D s £> 

• • 

OJ 

• 

v£> 

X 

X 

X 

X 

• 

X 

• 

X 

• 

X 

X 

X 

X 

X 

X 

X 

r- 

X 

• 

X 

• 

-O 

ft 

x 

M 

X 

«* 

V 

X 

e r 

pj 

•t 

'O 

J 

O' 

oO 

CD 

x> 

UJ 

p 

X 

ft 

X 

O' 

X 

n 

O' 

O'  P 

X 

X 

X)  p 

X 

r-4 

■X 

VP 

p 

T> 

P 

O' 

p 

X 

X 

X 

X 

X 

X 

X 

X 

X 

o 

X 

X 

l\J 

— • UJ 

O' 

VP 

<30 

UJ 

VP 

cr 

cr 

X 

X 

u 

X 

p 

CD 

ft 

X 

p 

p 

r> 

U' 

X 

-> 

u. 

O. 

JO  UJ 

X 

J 

O' 

X 

CJ 

•-* 

p 

O' 

UJ 

p 

X 

P 

U. 

f) 

X 

p 

L-J 

p 

ft) 

O' 

o 

CO 

• 

• • 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

X 

ft 

CD  © 

cu 

-/ 

JO 

o 

<p 

J 

«_ 

X 

p 

X 

X3 

ft 

ft 

X 

ft 

P 

p 

O' 

UJ 

X 

r> 

O 

X 

O 

IP 

-r 

-r 

X) 

r 

X 

X 

X 

X 

O 

ft 

O' 

cr 

P 

-J 

P 

X 

■> 

u. 

x 

J ft 

r 

•J 

-T 

JT 

X 

X 

J 

-J 

J 

T 

-r 

r? 

ft 

ft 

ft 

ft 

X 

p - 

f. . 

< 

u. 

Q. 

x 

<1 

• 

ft  r 

i_> 

»u 

X 

X 

l\J 

— • 

X. 

p 

IU 

u 

JJ 

X 

ft 

x> 

X 

cr 

O' 

p 

O' 

P 

X 

• 

© 

O'  X 

ip 

X 

-» 

a 

u. 

X 

X 

VP 

X 

© 

X 

»v 

© 

X 

X 

X 

X 

c 

X 

u 

►— 

f\j 

*-<  UJ 

C\J 

u> 

p- 

O' 

UJ 

O' 

VP 

<30 

X 

LD 

X 

X 

ft 

p 

X 

X 

p 

p 

P 

X 

u. 

ft 

x t 

IP 

X 

rr 

X 

« 

X 

O' 

VP 

P 

J 

X 

P 

P 

X 

CD 

o 

ft 

ft 

ft 

VP 

P> 

<J 

X 

P 

•ft  cr 

<r 

r 

vO 

PJ 

(V 

x 

X 

r 

X 

U 

u > 

X 

r 

X 

P 

X 

© 

p 

* 

r*’' 

UJ 

ft 

O 

«-< 

r **o 

fD 

•-< 

X 

IP 

o 

T 

X 

X 

X 

X 

p 

«- 

o 

X 

ft 

p 

ft 

p 

T 

UJ 

c 

ft  c 

a> 

C\J 

a 

H 

a 

c 

X 

UJ 

X 

p 

p 

"J 

X 

X 

X 

o 

X 

O 

P 

Z 

ft 

X 

o 

ip 

X f\J 

T) 

IP 

<c 

r 

X 

«J 

r 

p 

X 

X 

© 

P 

ft 

X 

© 

ft 

X 

X) 

x 

o 

o 

o 

X 

• • 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

»-4 

UJ 

— 1 

■D 

-X 

» J 

X 

X 

X 

X 

X 

X 

X 

X 

x 

X 

X 

X 

o 

© 

O' 

ft 

X 

ft 

ft 

ft 

X 

a 

u. 

_J 

x 

Pj  f\J 

PJ 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

*"• 

X 

X 

X 

t— 

o 

LT, 

o 

• 

ft 

L* 

LTV 

X 

P 

<M  -f 

a 

X 

vT 

X 

»-* 

VP- 

— < 

w* 

P 

© 

ft 

K 

X 

X 

CD 

X 

X 

r 

X 

_J 

P 

«t 

ft 

ip  cr 

<r 

O' 

or 

VP 

X 

cr 

X 

X 

p 

cr 

X 

X 

X 

© 

P 

X 

ft 

P 

cr 

o 

Or 

O 

IT 

X P 

lp 

X 

© 

X 

X 

— * 

X 

X 

X 

IT 

ft 

cr 

o 

p 

X 

ft 

o 

X 

p 

r> 

U 

u 

or 

iP  r\j 

w 

vp 

a 

w-4 

UJ 

O 

o 

f 

CD 

X 

ft 

t 

© 

© 

ft' 

ft 

■j 

-J 

X 

to 

P 

fV  x 

► 

X 

r»- 

X 

X 

r 

X 

u 

ip 

P 

P 

ir 

p 

X 

ft 

X 

CD 

© 

o 

o 

X 

CD 

sT 

y£> 

vP> 

vO 

-O 

sp 

U.' 

w 

u: 

w 

P 

u 

P 

P 

p 

p 

p 

X 

p 

P 

u 

_J 

• 

• • 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

<3 

u. 

-j 

X 

UJ 

n 

<r 

X 

• 

x 

O IP 

CD 

J 

*-< 

a 

cr 

X> 

UJ 

O' 

X 

ft 

X 

X 

X 

X 

X 

X 

O 

ft 

_j 

X 

« 

x 

-J  lO 

— t 

J* 

J 

CD 

© 

VT 

<p 

f 

O' 

X 

p 

X 

X 

X 

X 

ft 

CT 

CD 

ft 

a 

c 

M 

x 

r- 

U wP 

L 

U 

X 

X 

O'  cr 

O' 

© 

X 

X 

p 

ft 

u 

p 

P 

ft 

ft 

o 

ft 

UJ 

to 

o 

C-  o 

cr 

CD 

CD 

© 

CD 

o 

© 

*-4 

X 

X 

X 

X 

X 

X 

•-< 

X 

x 

X 

X 

r 

x. 

X* 

z 

X 

IT 

IP  IP 

ip 

O' 

u 

*P 

p 

p 

P 

tP 

U 

P 

p 

S 

p 

p 

ft 

a 

P 

u 

p 

U> 

u. 

UJ 

cr 

© 

CD  cr 

"3 

o 

r-> 

C- 

CD 

o 

o 

o 

o 

© 

CD 

© 

CD 

o 

© 

© 

o 

© 

© 

© 

O 

o 

•J 

• 

* 

x 

C.* 

X 

T 

X 

i 

<p 

o 

ft 

UJ 

IP 

CVJ  J’ 

C 

X 

CD 

p 

O' 

-r 

vP 

X 

C' 

ft 

K 

-» 

-T 

X 

ft 

ft 

© 

ft 

X 

© 

►-« 

cr 

x 

<r- 

X 

x> 

C 

-T 

X 

O' 

or 

X 

P 

cr 

X 

ft 

a 

p 

x 

p 

X 

X 

cr 

in 

o 

X 

ft 

«p  J 

<r 

o 

vO 

X 

X 

IT 

VP 

«D 

X 

p 

p 

© 

X 

ft 

ft 

© 

© 

O' 

O' 

u 

© O 

X 

x rvj 

cr 

•D 

-T 

cr 

s 0 

IP 

O' 

X 

X 

IP 

X 

K 

X 

-T 

r j 

p 

ft 

ft 

p 

u 

x 

z 

X 

PJ  x 

o 

t_> 

X 

O 

X> 

X 

ct 

O 

X 

P 

X 

ft 

X 

O 

O'  r 

X 

X 

r 

ft 

ft 

• 

»-• 

«r  sT 

J 

X 

© 

u 

u 

UJ 

X 

X 

O 

p 

X 

IP 

© 

O 

X 

p 

ft 

X 

X 

vO* 

r* 

i 

*3 

X 

r\>  r. 

X 

x 

t 

© 

O' 

O' 

or 

© 

X 

p 

P 

U' 

P 

► 

X 

CD 

x 

c 

c 

u 

UJ 

z 

a 

t-i 

X 

J ^ 

J 

J 

T 

X 

m 

r> 

X 

X 

X 

X 

X 

X 

ft 

ft. 

X 

ft 

ft 

u» 

n 

o 

<T 

> 

• 

• • 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

cr 

o 

S' 

X 

r 

*T 

V 

© 

p 

_> 

_J 

• 

cr 

o 

i 

r\> 

P*  Cl 

vx 

» 

u. 

UJ 

X 

X 

O IP 

X 

cr 

X 

X 

u> 

vP 

o 

IT 

X 

X 

X 

p 

p 

X 

ft 

x 

© 

p 

cr 

ft 

CD 

g 

cr 

»-4 

U 

VP  X 

U3 

LP 

O' 

«o 

VP 

JO 

u 

X 

x 

ft 

P 

p 

VP 

X 

J 

P 

CD 

X 

p 

M 

C. 

x 

X 

rr 

x x 

CJ 

X 

p 

p 

vf 

X 

vP 

© 

cr 

X 

CD 

X 

-» 

X 

p 

X 

p 

ft 

© 

r 

© 

Z 

3 

UJ 

U. 

u ^ 

X 

c 

■J 

<r 

P 

ip 

t5 

X 

X 

vr 

*“• 

X 

X 

p 

cr 

p 

o 

X 

X 

a 

p 

o 

u 

Cl 

Z 

ft 

x c. 

c 

X 

X 

UJ 

p, 

X 

r 

X 

* 

a 

r 

ft 

p 

4 

o 

p 

p 

ft 

ft 

UJ 

U J 

J 

U.  r 

cr 

s 

X 

X 

O' 

o 

© 

o 

© 

P 

X 

CT 

CT 

p 

p 

p 

© 

• 

a. 

X 

O 

X K) 

X 

X 

X 

X 

X x 

X 

X 

X 

X 

X 

X 

CD  ft 

X 

p 

© 

ft 

X 

P 

p 

*-u 

X 

u 

«3 

X 

X X 

X 

f 

X 

X 

X 

X 

r- 

X 

X 

X 

X 

X 

P 

P 

p 

p 

p 

p 

P 

X 

p 

« 

X 

«3 

-J 

• 

• • 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

o 

X 

O 

t 

cr 

a 

u 

«C 

u 

IT 

ft 

X 

o 

z 

to 

UJ 

x 

*1 

UJ 

X 

X 

-f  x 

•-< 

X 

© 

*-« 

ar 

X* 

vP 

O' 

© 

O' 

p 

« 

ft 

ft 

ft 

X 

ft 

ir 

X 

V 

JC 

ft 

c. 

-t  X 

X 

y£> 

J- 

X 

X 

cr 

-r 

cr 

X, 

X 

ft 

ft 

ft 

p 

ft 

ft 

o 

cr 

ft 

o 

C2 

• 

• • 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

* 

• 

• 

• 

• 

• 

• 

• 

e1 

u- 

ft 

J J 

vp 

VP 

vp 

vP 

vO  X 

X 

N 

o 

© 

ft  ft 

O' 

ft 

O' 

ft 

ft 

ft 

X 

X 

V 

cr.  u- 

IP  IP 

ip 

VP 

IP  IP 

IP  IP 

IP 

IP 

ip 

p 

IT 

P 

IX 

P 

VP 

p 

p »r 

P 

CJ 

X 

_J 

x 

u 

p 

CD  s£> 

X 

cr 

IP 

X 

X 

o 

X 

•D 

cr 

P 

X 

ft 

ft 

X X 

ft 

o 

X 

x 

H 

€_*  X 

X 

CD 

o 

vp 

o 

X 

P 

ft 

X 

X 

CD 

J 

ft 

X 

P 

O' 

© 

to 

-J 

«3 

• 

X ~4 

IP 

CD 

r^  X 

X 

vP 

O 

-* 

ft 

ft 

ft 

X 

p 

*4 

P 

O' 

P 

X 

X 

UJ 

• 

ft 

IP  wC 

sT 

X 

X 

X 

«> 

r 

cr 

ft 

© 

o 

© 

X 

X 

X 

X 

X 

X 

X 

ac 

f- 

o 

PJ  PJ 

X 

X 

X 

X 

X X 

X X 

X 

X 

X 

X 

ft 

ft 

ft 

X 

ft 

ft 

X x 

O 

PROBE  POSITION  (IN  FROM  CENTERBODY J 


■ 


Blade  Angle  65.0° 
Station  1 


PROBE  POSITION  UN  F RON  CENTERBODY) 


Blade  Angle  60.0' 
Station  1 


73 


I 


I 


PROBE  POSITION  IN  PROM  CENTERBODY) 


Blade  Angle  60.0° 
Station  1 


Figure  5?  Flow  VS  Probe  Position 
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TOTAL  FLCM  RATF  IS  13DlO<.6<.  IVSEC 


PROBE  POSITION  UN  PROM  CENTERBODY) 


Blade  Angle  69.7° 
Station  3 


54-0  56-0  58-0  60-0  62-0  64-0  66-0 

FLOW  ANGLE  (DEG) 

Figure  55  Flow  Angle  VS  Probe  Position 
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Blade  Rngle  65.0° 
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PROBE  POSITION  (IN  FROM  CENTERBODY) 


Blade  Angle  60.0° 
Station  3 


FLOW  ANGLE  (DEG) 


Figure  63  Flow  Angle  VS  Probe  Position 
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Figure  64  Flow  VS  Probe  Position 
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Calculated  results  are  then  shown  in  Tables  6 through  11  with  correspond- 
ing plots  of  mass  flow  distribution  immediately  following  the  tables  in 
Figures  44,  48,  52,  56,  60,  and  64.  The  division  of  the  annulus  area 
into  circumferential  strips  for  mass  flow  calculation  is  shown  in 
Figure  65.  A computer  program  was  written  to  generate  these  answers  and 
complete  details  of  the  calculations  are  shown  in  Appendix  J.  A total  of 
one  hundred  eleven  (111)  different  traverses  were  made  during  the  course 
of  this  experimental  program  and  the  complete  tabulation  of  raw  data  was 
transferred  to  computer  cards  at  the  Air  Force  Aero  Propulsion  Laboratory. 
A thorough  review  of  all  these  runs  and  a series  of  repeat  runs  produced 
the  conclusion  that,  for  any  given  cascade,  the  results  were  both  consis- 
tent and  repeatable.  Therefore,  those  shown  here  are  deemed  representa- 
tive and  adequate  to  establish  the  correct  conclusions.  I r,  all  cases,  due 
to  the  location  of  pressure  taps  on  the  probe  and  the  need  for  perating 
clearance  between  the  centerbody  and  the  end  of  the  probe,  the  innermost 
reading  was  approximately  0.3  inches  above  the  surface  of  the  centerbody. 
At  the  outside  wall,  due  to  the  liner  penetration  scheme,  the  probe  could 
be  withdrawn  such  that  the  pressure  taps  were  physically  at  the  wall. 

Thus,  in  all  plots,  the  boundary  layer  effect  at  the  outer  wall  is  clearly 
visible  while  at  the  inner  wall  the  thickness  of  the  boundary  layer  is 
generally  less  than  the  0.3  inch  limit.  The  fact  that  near  the  outer  wall 
there  is  a stronger  boundary  layer  effect  than  near  the  hub,  can  be 
attributed  to  the  flow  path  following  the  bellmouth  entrance.  The  accumu- 
lated boundary  layer  growth  along  the  entrance  bellmouth  is  manifested  in 
the  greater  thickness  of  the  boundary  layer.  Obviously,  the  same  effect 
took  place  near  the  centerbody  but  on  a much  weaker  scale.  So  each  plot 
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MASS  FLOW 
AREA  DIVISION 


can  be  considered  as  demonstrating  several  zones;  i.e.,  outer  wall 
boundary  layer,  fully  developed  channel  flow,  distorted  flow  due  to  the 
effects  of  increasing  blade  blockage  at  inner  radii  and  finally  inner 
wall  boundary  layer  effect  which  could  not  be  measured  but  is  certainly 
there.  As  can  be  seen  readily  from  the  total  pressure  and  mass  flow 
distribution  plots,  uniformity  across  the  channel  has  been  achieved. 

The  scales  selected  for  these  plots  tend  to  exaggerate  what  in  reality 
are  relatively  small  percentage  changes.  For  example,  taking  the  worst 
case  shown  in  Figure  41,  there  is  a maximum  variation  of  approximately 
368  to  398  inches  of  water  absolute  total  pressure.  Therefore,  the  maxi- 
mum change  of  30  inches  compared  to  an  average  value  of  388  inches  gives 
a percentage  variation  of  + 3.87%.  These  scales  were  consciously 
selected  to  demonstrate  both  the  sensitivity  of  the  instrumentation  used 
and  the  repeatability  of  the  data.  A sample  typical  repeat  run  comparison 
is  shown  in  Figures  66,  67,  68,  and  69,  for  the  65°  blade  row.  In  this 
case,  the  maximum  deviation  between  runs  of  total  pressure  is  approximately 
5 inches  of  water  out  of  a nominal  absolute  value  of  395.  Thus,  the 
percentage  change  is  approximately  +_  0.633%.  The  flow  distribution  com- 
parison shown  in  Figure  69  indicates  excellent  repeatability  and  the 
legend  gives  the  receiver  pressures  for  the  respective  runs.  These  were 
run  on  different  days  and  the  receiver  pressure  was  always  a function  of 
facility  controls  and  the  atmospheric  pressure. 

This  uniformity  of  total  pressure  and  flow  per  unit  area  of  the 
channel  is  the  essence  of  the  proof  of  Engelman's  analysis.  Achievement 
of  uniform  axial  mass  flow  distribution  at  high  subsonic  Mach  numbers 
with  nontapered,  nontwisted  blades  is  not  predicted  by  the  free  vortex 
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Blade  Angle  65.0° 
Station  1 
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Blade  Angle  65.0° 
Station  1 


Figure  68  Flow  Angle  VS  Frobe  Position 
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theory.  The  compressibility  effects  are  clearly  real  and  the  analysis 
by  Engelman  is  better  able  to  predict  actual  flow  parameters.  Comparative 
plots  are  presented  later  in  this  chapter  to  illuminate  this  point  in 
Figures  109  and  110. 

As  an  aid  in  relating  the  data  taken  to  the  geometry  of  the  flow 
path  in  the  cascade, several  graphical  presentations  and  plotting  schemes 
were  exercised.  These  were:  (1)  velocity  vectors  with  axial  and 
tangential  components  across  representative  traverses,  (2)  three  dimen- 
sional plots  of  velocity  profiles  across  the  circumferential  range  inves- 
tigated, and  (3)  plots  of  constant,  as-measured  flow  angle  across  the 
annulus.  These  have  been  collected  together  in  groups  for  ease  in 
reference  and  are  shown  in  Figures  70  through  108.  Tables  12  and  13 
organize  the  collection  of  velocity  vectors  and  components  presented. 
Discussion  of  these  presentations  follows  each  grouoing  as  appropriate. 


10°  Right  of  T.E. 
65°  Blades 
Station  1 


Figure  71  Velocity  Vectors  & Components 


2.5°  Right  of  T.F. 


60°  Blades 
Station  1 


Figure  72  Velocity  Vectors  & Components 


7.5°  Right  of  T.E. 


65°  Blades 
Station  3 


Figure  74  Velocity  Vectors  & Components 
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2.5°  Left  Figure  95 

4.5°  Left  Figure  96 


65°  Blades 
Station  1 


Figure  83  Velocity  Vectors  & Components 


2.5°  Right  of  T.E. 
65°  Blades 
Station  1 


Figure  84  Velocity  Vectors  S Components 


Station  3 


Figure  89  Velocity  Vectors  4 Components 


Station  3 


Figure  90  Velocity  Vectors  & Components 


Station  3 


Figure  91  Velocity  Vectors  & Components 


2.5°  Left  of  T.E.  ^ 0,o 

60°  Blades 
Station  3 
Low  Flow 


Figure  95  Velocity  Vectors  & Conponents 


4.5°  Left  of  T.E.  ' 0 

60°  Blades 
Station  3 
Low  Flow 

Figure  96  Velocity  Vectors  & Components 


2 


The  velocity  vectors  and  their  axial  and  tangential  components  plotted 
in  Figures  70  through  96  include  the  entire  traverse  for  each  case. 

It  should  be  noted  that  in  all  these  plots,  the  radial  location  is  shown 
in  inches  from  the  test  assembly  centerline  rather  than  from  the  center- 
body  surface.  Once  again,  the  outer  wall  boundary  layer  effect  is 
clearly  shown  while  the  inner  wall  or  centerbody  layer  generally  does 
not  even  start  to  show  up  due  to  the  last  reading  being  approximately 
0.3  inch  above  the  surface.  With  the  various  axes  labeled  and  graduated, 
an  absolute  feel  for  relative  velocities  at  each  radial  station  can  be 
easily  obtained.  In  addition  it  is  possible,  with  a little  more  effort 
at  visualization,  to  ascertain  the  relative  changes  in  flow  angle  at 
each  radial  station. 
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Station  1 - 21  Traverses 
69.7°  Blade  Anale 
High  Flow 


Fiaure  97  3D  Composite  of  Velocity  Profiles 


Station  3-8  Traverses 
69.7°  Blade  Angle 
High  Flow 


Figure  98  3D  Composite  of  Velocity  Profiles 


Station 

65 


1 - 16  Traverses 
Blade  Anqle 
High  Flow 


Figure  99  3D  Composite  of  Velocity  Profiles 


Station  3-7  Traverses 
60c  Blade  Angle 
Low  Flow 


Figure  101  3D  Composite  of  Velocity  Profiles 


The  three-dimensional  composites  of  velocity  profiles  plotted  in 
Figures  97  through  101  incorporate  representative  runs  in  each  case  over 
the  entire  range  of  channel  traverses  recorded  in  this  experiment.  In 
order  to  prevent  complete  obliteration  of  the  views  by  including  all 
velocity  vectors,  only  the  hub  and  tip  vectors  were  printed.  Then  the 
end  points  of  all  the  others  were  printed  to  generate  the  profiles  shown. 
The  only  purpose  of  these  plots  was  to  attempt  to  aid  visualization  of  the 
flow  field  actually  measured  and  they  are  not  intended  to  be  scalable 
directly;  i.e.  the  profiles  are  correctly  plotted  with  respect  to  each 
other  but  on  an  absolute  velocity  basis  cannot  be  scaled. 

When  comparing  the  station-one  plots  with  the  station-three  plots, 
it  becomes  apparent  that,  as  expected,  the  blade  wake  effects  are  start- 
ing to  mix  out  and  while  still  there,  are  much  reduced.  However,  any 
assumption  of  truly  wake  free  flow  is  still  not  valid  even  at  station 
three.  Due  to  the  limited  circumferential  travel  available  in  the  test 
assembly,  the  actual  wake  flow  could  not  be  tracked  at  station  three 
from  the  specific  blades  surveyed  at  station  one.  The  wake  effects  seen 
here  at  station  three  are  generated  by  other  blades,  and  in  order  to  make 
any  claim  to  absolute  tracking  of  wake  decay,  one  must  assume  that  iden- 
tical blade  profile,  surface  roughness,  blade  angle,  and  blade  spacing 
exist  for  every  blade  in  the  cascade.  Due  to  the  intense  quality  control 
exercised  during  the  fabrication  of  the  test  assembly,  this  is  considered 
a reasonable  approach  with  errors  induced  from  these  sources  to  be  less 
than  any  flow  measurement  errors. 
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Figure  102  60°  Blades,  Sta  1 

Constant  Flow  Angle  Lines 


Figure  103  60°  Blades,  Sta  3,  18.10"  Hq 
Constant  Flow  Angle  Lines 


Lines  of  constant  flow  angle  across  the  annulus  are  plotted  in 


Figures  102  through  108.  These  plots  make  use  of  all  data  taken  and  the 
circumferential  range  surveyed  in  each  case  is  clear.  For  reference 
to  the  blade  row,  the  tick  marks  at  the  hub  and  tip  lines  represent 
blade  trailing  edges.  The  pressure  readings  shown  in  the  titles  of 
Figures  103, 104,  andl05are  the  receiver  pressures  and  indicate  less 
than  maximum  flow  conditions;  i.e.  what,  in  general,  were  called  "low  flow" 
runs  during  the  course  of  this  investigation.  These  runs  were  made  to 
gather  comparative  data  to  test  the  notion  discussed  throughout  this 
dissertation  relative  to  the  free  vortex  being  most  valid  at  lower  Mach 
numbers  while  the  Engelman  analysis  would  be  most  valid  at  high  Mach 
numbers.  In  general,  the  "high  flow"  runs  were  conducted  with  a receiver 
pressure  of  11  inches  of  Hq  absolute  or  less. 

In  these  plots,  the  general  blade  effects  can  be  seen  as  well  as 
the  variation  in  effect  across  the  channel.  The  lines  are  labeled  in 
degrees  and,  to  prevent  clutter,  several  interior  islands  and  lines  were 
not  specifically  marked.  However,  simply  note  that  all  lines  were 
indexed  one  degree  apart  so  that  any  unmarked  line  can  be  identified  by 
referring  to  a nearby  marked  one.  One  word  of  caution  in  interpreting 
the  plots  is  in  order;  i.e.  due  to  the  nature  of  the  plotting  program, 
if  no  data  exist  in  an  interval,  it  merely  draws  a straight  line 
between  the  points  available.  Therefore,  these  effects  must  be  evaluated 
accordinoly  when  using  the  plots,  but  the  presentation  in  this  fashion 
was  considered  worthwhile  as  still  anotner  wav  to  aain  insight. 

A final  observation  which  can  be  made  by  comparing  the  plots  for 
station  three  with  those  for  station  one  is  that  there  is  very  little 
reduction  in  flow  angle.  The  axial  distance  between  stations  one  and 
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three  is  four  inches  or  approximately  3.2  chords. 

In  particular,  when  comparing  the  69.7°  data  shown  in  Figures  107  and 
108,  the  central  channel  flow  shows  no  real  change  at  all,  and  the  outer 
wall  boundary  layer  appears  to  be  growing  very  slowly. 

This  growth  is  evidenced  by  the  reduction  in  radius  at  which  the  flow 
angle  starts  to  be  reduced  by  the  wall  boundary  layer  effect. 

In  summary,  by  referring  to  one  or  more  of  these  graphical  presenta- 
tion schemes,  it  is  possible  to  get  a clear  visualization  of  the  measured 
flowfield,  including  a reasonable  feel  for  the  absolute  value  of  the  flow 
angles  and  velocities. 

The  Reynolds  numbers  were  calculated  for  representative  high  and  low 
flow  runs  over  the  entire  flow  range  investigated  with  each  set  of  blades 
and  these  results  recorded  in  Appendix  I.  In  each  case,  both  the  highest 
and  lowest  number  achieved  are  shown.  The  flow  ranges  investigated  were 
as  shown  in  Table  14  and  the  corresponding  Reynolds  numbers  are  included 
for  reference.  The  characteristic  dimension  used  to  calculate  these 
Reynolds  numbers  was  twice  the  channel  height  (blade  span). 


Table  14 

Flow  Ranges  Investigated 


Blade 

Angle 

Flow  Rate  Range 

lb/sec 

Reynolds  No  Range 

( io'6) 

69.7° 

11.6  - 14.4 

1 .416  - 1 .957 

65° 

12.8  - 15.5 

1.312  - 1.911 

60° 

13.5  - 16.0 

1.151  - 1.566 

These  rates  were  calculated  based  on  the  traverse  data  collected,  and 
typical  sets  of  results  are  shown  in  Tables  6 through  11.  The  annulus 
division  into  strips  for  these  calculations  is  shown  in  Figure  65. 

Surface  roughness  measurements  for  all  parts  of  the  flow  path  were 
made  with  a portable  profilometer  and  these  results  are  recorded  in 
Appendix  H.  The  flow  path  boundaries  are  machined  aluminum  (annulus 
inner  wall),  epoxy-varnished  mahogany  (annulus  outer  wall)  and  extruded 
aluminum  (blades).  In  the  case  of  the  blades,  the  extrusion  die  marks 
were,  of  course,  perpendicular  to  the  flow.  In  qeneral , however,  all  parts 

of  the  flow  path  were  quite  smooth  mechanically  and  based  on  the  work  of 

( 95) 

Bammert  and  Sandstede'  , can  be  considered  hydraulically  smooth  at  the 
Reynolds  numbers  achieved  in  this  experiment. 

Calculations  were  then  made  to  compare  these  measured  results  with 
various  theoretical  treatments  and  the  details  of  that  phase  are  covered 
in  Chapter  IV. 
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CHAPTER  IV 


COMPARISON  WITH  THEORETICAL  CALCULATIONS 

Theoretical  flow  angle  calculations  were  made  in  the  following  four 
ways  to  generate  predictions  which  could  be  compared  with  the  measured 
results  from  this  experiment: 

(1)  Free  vortex  assumption 

(2)  Engelman  analysis 

(3)  Two  dimensional  computer  solution 

(4)  Three  dimensional  computer  solution 

All  four  schemes  are  discussed  in  detail  below  and  the  comparisons  are 
made  following  the  calculation  descriptions.  The  comparative  results  are 
presented  in  graphical  and  tabular  form  as  appropriate. 

First,  a computer  program  was  written  to  find  the  theoretical  flow 
angles  as  a function  of  pressure  ratio  and  as  predicted  by  both  the 
free  vortex  analysis  and  the  Engelman  analysis.  The  angle  : * is  defined 
as  the  flow  angle  where  Mach  number  equals  one. 

The  governing  equation  for  angles  as  predicted  by  Engelman's 
analysis  is: 

_ k+1 

c°s:+  , y 
cos : 


k+1 


2 + (k-l)M 


2TTTF 


(4.1) 
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This  equation  can  be  obtained  by  use  of  the  following  method: 

Substituting  the  equation  for  the  geometric  relationship  as  defined  by 
Figure  A-l  in  Appendix  A 


Vz  = Vcos? 

into  the  continuity  equation 


and  rearranging 

cos  ;♦  _ V 

COS  9 ~ 

Substituting  the  following  isentropic  relationships 


(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 


(4.7) 


(4.8) 


£_  = £ 1 + ( Jill  ) M2  J - p (4.9) 

into  equation  (4.7)  and  evaluating  at  M*  = 1 , equation  (4.1)  is  obtained. 
In  this  equation  M is  the  Mach  number  and  t is  the  predicted  angle. 


143 


Then,  assuming  k = 1.4,  and  solving  equation  (4.1)  for  4 and  equation 
(4.9)  for  M gives 


4 = cos 


M = 4/5 


, r cos  (2  + 0.4  m2\  3i 

L M \ 24  / J 


Finally,  combining  (4.10)  and  (4.11)  gives 


cos  4* 


4 = COS 


I'M 

Lw 


■ • • »!(;) 


(4.10) 


(4.11) 


(4.12) 


FREE  VORTEX  ANALYSIS 

For  the  Free  Vortex  Analysis,  the  governing  equations  for  flow  angles 
obtained  in  Appendix  A are: 


sin 


(4.13) 


where. 


(4.14) 


(k+1)  P'  1 
o 

(k-l)sin  4^ 


Combining  (4.131  with  (4.14)  for  k = 1.4  gives: 


<■  = sin 


/ 6p • ^857_5  * 

sin 

it/1  - • 2 

/ c /•,  d0.2857\ 

M sin  4* 

v ( ) 

( 

L"  \ / J 

« 

(4.15) 


Equations  (4.12)  and  (4.15)  are  both  functions  of  P and  In  Figures 

109  and  110,  the  flow  angles  are  plotted  as  a function  of  channel  posi- 
tion. A value  of  was  calculated  and  used  as  a constant  parameter. 

Then  the  measured  pressure  ratio  at  a particular  channel  position  was 
used  to  find  the  theoretical  { values  for  both  the  Free  Vortex  analysis 
and  the  Engelman  analysis.  The  calculation  of  for  Figure  109  was  made 
as  follows: 

continuity  of  axial  mass  flow  requires, 

cVz 

so  that 


*V*  cos  :* 
V, 


cos 


- 


p*V< 


(A. 32) 

(A. 16) 


Using  the  mid-channel  data  obtained  at  this  station  and  the  isentropic  com- 
pressible flow  functions  from  Table  30  of  the  Gas  Tables  by  Keenan  and 
Kay^37,  pp'  139,  a value  for  cos  can  be  found  as  follows: 

At 


M = 1 — * 0.63394 


and  at 


V = 0.7 


= 0.79159 


so 


£_  = • 791 59  _ i 249 

.63394  ,c 
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then 


1.249  V. 


= cos  <?* 


From  the  reduced  data. 


Vz  = 362  ft/sec  and  Tt  = 542°R 


Again,  from  the  Gas  Tables, 


so 


Hence 


at  M = 0.7  j = 0.91075 
't 


T = (0.91075)  (542)  = 493. 6°R 


V*  = 49.02  T = 493.6  (ref)  = 1089  ft/sec 


Substitution  into  eq  (4.16)  yields 


thus 


cos 


(1.249)  (362) 
' — '1089 


0.41519 


A 

♦ ★ 


65.5° 


A similar  calculation  for  conditions  of  Figure  (110)  yields  t*  = 70°. 

As  can  be  seen  from  these  plots,  the  Engelman  solution  more  nearly 
approximates  the  measured  real  flow  field  at  high  Mach  numbers  than 
does  the  free  vortex  solution. 
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PROBE  POSITION  (INCHES  FROM  CENTERBODY) 


Next,  in  order  to  permit  a comparison  of  the  data  resulting  from 
tnis  experiment  with  current  computerized  analytical  prediction  techniques, 
both  two-dimensional  and  three-dimensional  treatments  were  applied  to  the 
cascade  flow  path.  These  programs  were  developed  under  contract  to  the 
Air  Force  Aero  Propulsion  Laboratory  Turbine  Engine  Division  and  are  being 
routinely  applied  to  advanced  design  efforts.  As  with  all  analytical 
schemes,  certain  assumptions  must  be  made  to  permit  the  calculations  to 
proceed  or  to  minimize  computer  time  requirements.  Therefore,  the  appli- 
cation of  these  programs  to  this  experimental  work  was  an  excellent  addi- 
tional opportunity  to  compare  predictions  with  the  measured  data.  The 
two-dimensional  program  was  prepared  by  Wysong  1 of  the  General  Electric 
Company  and  the  three-dimensional  program  was  prepared  by  Dodge' at  the 
Garrett  Corporation.  The  results  of  this  activity  are  presented 
sequentially  below. 

The  Turbine  Design  System  (TDS)  developed  by  Wysong'  at  the 
General  Electric  Company  includes  various  modules,  each  of  which  specifi- 
cally addresses  a major  component  of  the  overall  design  process.  Its 
primary  overall  purpose  is  to  provide  a preliminary  design  tool  and  is 
set  up  as  an  interactive  scheme  intended  for  use  on  a terminal  equipped 
with  a cathode  ray  tube  and  a hard-copy  device.  For  the  purpose  of  com- 
parison with  the  data  from  this  experiment  the  CASC  Module  was  exercised 

for  two  representative  cases.  The  following  description  of  the  workings 

(26) 

of  this  module  is  quoted  directly  from  Wysong'  : 

"The  CASC  module  performs  a blade-to-blade  flowfield  analysis 
on  a cascade  of  airfoil  contours  for  the  purpose  of  determining 
Mach  number  and  pressure  distributions  on  the  airfoil  surfaces. 
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These  distributions  provide  the  basic  ingredients  for  evaluating 
blade  section  performance  and  for  predicting  blade  heat  flux 
rates. 

"The  streamline  curvature  technique  (SLC)  is  used  to  solve  the 
steady,  compressible,  inviscid  flow  equations.  The  calculation 
domain  (Fig.  8)  is  a lamina  of  varying  thickness  conforming  to 
the  shape  of  an  axisimnietric  meridional  stream  surface  which 
cuts  the  blade.  Within  this  lamina,  a number  of  confluent 
stream-tubes  are  bound  together  to  obtain  the  total  flow  in  the 
blade-to-blade  passage.  This  binding  of  stream-tubes  is  accom- 
plished by  a grid  of  streamlines  and  orthogonal  lines  which 
provides  the  basic  structure  for  the  computing  procedure.  SLC 
has  been  found  to  be  an  accurate  method  even  for  coarse  grids. 
The  accuracy  accrues  from  the  SLC  formulation,  which  is  imbedded 
in  a natural  intrinsic  coordinate  system,  and  is  enhanced  by 
the  fact  that  the  cross-stream  velocity  distribution  utilizes 
airfoil  surface  curvatures  as  boundary  information. 

"The  CASC  version  of  SLC  employs  "type-dependent  differencing" 
to  account  for  subsonic  and  supersonic  flow  regimes  and  a 
direct  matrix  solution  for  reliable  convergence.  Repeating 
boundary  conditions  are  enforced  outside  of  the  "covered" 
blade-to-blade  passage.  The  governing  equations  include  cen- 
trifugal and  Coriolis  force  terms,  hence,  both  stationary  and 
rotating  blade  rows  can  be  treated. 

"A  brief  description  of  the  solution  is  as  follows.  An  initial, 
approximate  grid  of  streamlines  and  orthogonals  is  established. 
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Streamline  curvature  is  evaluated  at  each  grid  point.  The 
Crocco  form  of  the  momentum  equation  is  integrated  along 
orthogonals  to  obtain  cross-stream  velocity  variations  and  the 
continuity  equation  is  integrated  along  orthogonals  to  determine 
corrected  streamline  positions  for  the  current  curvature  field. 
Entropy  and  rothalpy  are  held  constant  in  the  streamwise 
direction.  Grid  point  adjustments  are  calculated  by  considering 
the  difference  between  computed  and  current  streamline  positions 
and  the  effect  of  the  implied  curvature  modification  on  the 
integrated  momentum  equation.  The  streamlines  are  then  reposi- 
tioned according  to  the  calculated  adjustments.  The  process  is 
repeated  until  the  maximum  calculated  point  adjustment  is  less 
than  a prescribed  tolerance. 

"Since  the  movement  of  any  one  grid  point  alters  the  velocities 
at  nearby  points  through  a change  in  curvatures,  the  point 
adjustments  throughout  the  field  are  made  simultaneously  through 
a direct  matrix  solution.  This  departure  from  classical  stream- 
line curvature  methods  accounts  for  the  rapid  convergence 
characteri sties  of  SLC. 

"Input  to  CASC  falls  into  two  categories.  The  first  is  geo- 
metric data  which  include  grid  formation  parameters,  blade 
section  surface  coordinates  and  slopes,  and  meridional  surface 
definition,  all  of  which  come  from  the  BLDDES  module  results. 

The  second  category  is  aerothermodynamic  data  which  include  gas 
angles,  inlet  Mach  number,  inlet  total  pressure,  inlet  total 
temperature,  and  downstream  static  pressure.  All  of  these  come 
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from  the  ATREC  module  results.  They  are  part  of  the  data 
extraction  function  of  BDPREP  and  are  passed  along  through  the 
BLDDES  file  for  availability  to  CASC. 

"Execution  of  the  module  involves  accessing  one  blade  section 
at  a time  from  the  input  file.  Terminal  queries  and  responses 
permit  alteration  of  grid  formation  parameters  when  unusual 
circumstances  dictate  changes  to  the  stereotyped  data  supplied 
from  BLDDES.  Other  input  items  can  also  be  changed  to  guide 
the  course  of  the  iterations.  Output  at  the  terminal  includes 
plots  of  the  grid  structure  and  of  the  blade  surface  Mach  num- 
ber distributions  (Fig.  9).  Optional  printouts  of  diagnostics 
and  solution  results  a-e  available.  These  results,  coupled 
with  the  designer's  expertise,  are  used  to  judge  the  aerody- 
namic performance  of  the  blade  sectior  and,  if  deemed  necessary, 
to  determine  where  contour  adjustments  should  be  made  when 
rerunning  the  section  through  the  BLDDES  module. 

"This  procedure  continues  until  all  blade  sections  of  the  input 
file  have  been  evaluated.  Then  the  ELDDES  module  is  reactivated 
to  create  new  versions  of  disputed  sections  and  to  produce  a 
new  input  file  for  CASC.  CASC  is  rigged  to  read  the  file 
selectively,  so  only  the  altered  sections  are  treated.  The 
final  output  file  from  CASC  has  data  sets  containing  updated 
input  variables  and  calculation  results  for  at  least  one 
acceptable  version  of  hub,  pitch,  and  tip  design  sections. 
Multiple  versions  of  one  or  more  of  these  sections  may  be 
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Calculations  were  made  for  five  sections  along  the  span  of  the 
blades  at  the  following  radii:  5.5  (Root),  6.125,  6.75  (Pitch),  7.375, 
and  8.00  (Tip).  The  results  of  these  calculations  for  the  69.7°  blade 
row  (48  blades)  are  summarized  in  Table  15,  while  those  for  the  60c  blade 
row  (54  blades)  are  shown  in  Table  16  . 

Table  15 

Flow  Angle  Comparison 
(48  Blades)  69.7° 


61.09 

59.71 

59.95 

60.36 

58.90 
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Before  making  a comparison  with  the  data  collected  in  this  experiment 
recall  that: 

(1)  This  analysis  is  based  on  a two  dimensional,  inviscid  approach. 

(2)  The  closest  data  available  from  this  experiment  was  taken  0.75 
inches  (0.58  chord)  behind  the  exit  plane. 

Now,  the  quickest  way  to  compare  results  is  to  refer  to  Figure  107  for 
the  69.7°  blade  row  and  Figure  102  for  the  60c  blade  row.  Both  these 
figures  are  in  Chapter  III,  Results.  To  summarize  such  a comparison, 
it  appears  that  the  analysis  is  consistently  predicting  a flow  angle 
several  dearees  higher  than  that  measured.  Fnr  the  69.7°  blades, 
the  hub-tip  ranqe  measureo  was  58°-64c‘  vs  the  prediction  of  71.5°- 
69. 2C.  It  is  also  considered  significant  that  the  prediction  is  inverted 
compared  to  the  measurements;  i.e.,  the  angle  was  predicted  to  decrease 
from  hub  to  tip  while  the  measured  values  increased  from  hub  to  tip. 

For  the  60°  blades  the  hub-tip  range  measured  was  53c-57  vs  the 
prediction  of  61.09-58.90;  again,  the  predictions  were  inverted  compared 
to  the  measurements  anc  several  degrees  too  high  in  all  cases. 

A complete  description  of  the  theory  and  numerical  procedures,  used 
in  the  three-dimens ional  calculations,  is  included  in  the  AFAPL-TR-77-64 , 
"3-D  Heat  Transfer  Analysis  Program,"  October  1977  (Dodge^^).  This 
calculation  scheme  basically  combines  a relaxation  solution  technique 
with  a marching  solution  in  accordance  with  the  logic  diagram  shown  in 
Figure  1 1 1 taken  from  Dodae^’  To  quote  Dodge^1’  ^ specif- 

ically on  the  algorithm: 

"Three  separate  numerical  algorithms  are  required.  They  are: 
(a)  the  setup  of  the  geometry  and  its  conversion  coefficients. 
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Figure  111  Program  Logic  Diaoram 
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(b)  a relaxation  solution  for  potential,  and 

(c)  a marching  solution  for  temperature  and  the  three  compo- 


nents of  viscous  velocity." 

The  results  of  applying  this  technique  to  the  69.7°  blade  row  (48 
blades)  are  shown  in  Figure  112  where  L and  N are  the  radial  and  circum- 
ferential coordinates  respectively.  As  can  be  seen  there,  the  computer- 
generated output  plot  gives  predicted  flow  angles  at  the  trailing  edge 
plane  of  the  blade  row.  Now,  as  with  the  two-dimensional  discussion, 
probably  the  quickest  way  to  compare  measured  results  to  these  predic- 
tions is  to  use  Figure  107  in  Chapter  III. 

In  this  case,  for  the  69.7°  blades  the  hub-tip  range  measured  was 
again  58°-64*  versus  the  three-dimensional  prediction  of  approxima tely 
68°-72°  on  the  pressure  side  of  Figurell2  and  a nearly  constant  66r  on 
the  suction  side.  The  predicted  midsoan  ranqe  across  the  blade  passage 
exit  ranges  from  70°  on  the  pressure  side  to  66°  on  the  suction  side. 

Tne  measured  mid-span  range  from  Fiqure  107  is  61c-63°.  It  is  recog- 
nized that  the  reader  can  view  and  interpret  the  Figure  112  results  at 
will,  based  on  his  experience,  to  arrive  at  different  bounds  on  these 
ranges  depending  on  which  zones  he  chooses  to  ignore. 
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Now,  to  rationalize  this  apparent  discrepancy  of  several  degrees 
between  predictions  and  measurements,  two  possible  explanations  are  con- 
sidered worthy  of  note  here: 

(1)  The  real  blade  exit  angles  are  higher  than  the  downstream 
measured  angles  at  station  one  and  that  in  0.75  inch  the  flow  "untwists" 
or  straightens  several  degrees;  i.e.,  the  predictions  are  in  fact  correct 
and  so  are  the  measurements. 

(2)  The  real  blade  exit  angles  do  not  change  measurably  at 
these  Mach  numbers  in  such  a short  distance;  i.e.,  the  predictions  are 
actually  in  error  by  the  several  degrees  indicated  by  a simple  comparison. 

To  believe  (1),  it  must  be  accepted  that  a mechanism  exists  to  cause 
an  instant  straightening  of  the  flow  upon  leaving  the  blade  row  which  then 
was  measured  at  station  one  and  remained  practically  constant  for  at  least 
another  4 inches  (3.2  chords)  downstream  to  be  measured  aaain  at  station 
three.  Again,  it  must  be  rememoered  that  the  predictions  are  for  the 
exact  trailing  edqe  plane  while  the  nearest  experimental  data  was  0.75 
inch  (0.58  chord)  downstream  from  the  trailing  edge.  The  exact  amount 
of  flow  straightening  that  occurs  at  these  Mach  numbers  in  that  axial 
distance  is  debatable.  In  addition,  it  has  already  been  pointed  out, 
while  discussing  the  constant  flow  angle  plots  in  Figures  102  through 
108,  that  little  change  in  measured  angle  occurs  between  station  one 
and  station  three.  Furthermore,  shown  in  Figure  K-3for  the  65c  blade 

row  is  a plot  of  flow  angles  measured  across  the  channel  during  a check 
run  at  station  7 which  is  another  eight  inches  (6.2chords)  downstream. 

As  can  be  seen  there,  the  hub-to-tip  range  of  angles  varies  from  approxi- 
mately 52°  to  58°.  Now,  referring  to  Figure  106  which  plots  constant 
flow  angle  lines  for  the  65°  blades  at  station  1,  it  can  be  seen  that  the 
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measured  hub-tip  range  is  approximately  54°-62°.  Therefore,  even  in  a 
distance  of  12.75  inches  (9.86  chords)  downstream  from  the  blade  exit 
plane  there  is  only  approximately  2 to  4 degrees  of  "untwist"  in  the  flow. 
This  comparison  clearly  lends  credence  to  the  notion  that  little  flow 
straightening  is  occurring  downstream  at  these  Mach  numbers.  For  reference, 
the  complete  set  of  data  taken  and  flow  calculations  made  for  this  check 
run  are  shown  in  Appendix  K. 

To  believe  (2),  it  must  be  accepted  that  the  exit  flow  angles 
produced  at  these  Mach  numbers  will  never  exceed  the  metal  angles  and  that 
no  effects  other  than  boundary  layer  buildup  exists  to  cause  straightening 
of  the  flow.  This  author  subscribed  to  (2),  and  both  the  following  con- 
clusions about  the  three-dimensional  predictions  and  the  previous  dis- 
cussion of  two-dimensional  results  were  tempered  by  that  opinion. 

At  this  point,  strictly  on  the  basis  of  gross  flow  angle  Drediction, 
it  is  difficult  to  make  a case  that  the  three-dimensional  is  considerably 
better  than  the  two-dimensional.  If  the  two-dimensional  procedure  were 
applied  enough  times  across  the  span  of  the  blades,  one  could  expect  to 
generate  a still  better  picture  of  the  flow  albeit  still  inverted.  How- 
ever, it  would  appear  that  the  real  strenath  of  the  three-dimensional 
effort  lies  in  the  detail  available  with  regard  to  secondary  flows  and 
boundary  effects  which  the  two-dimensional  approach  will  never  yield  no 
matter  how  thin  the  spanwise  slices  become.  The  overall  result  pursued 
by  Dodge^',  after  all,  was  a heat  transfer  analysis  program  applicable 
tu  turbine  design.  Thus,  the  final  conclusion  regarding  utility  of  the 
three-dimensional  scheme  is  that  it  certainly  produces  more  detailed 
results  across  the  entire  blade  passage  and,  in  the  hands  of  an  experienced 
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designer  with  access  to  adequate  computer  time,  could  eliminate  the 
requirement  for  much  experimental  work.  However,  final  testing  and 
hardware  iteration  will  likely  always  be  required  if  maximum  performance 
is  the  goal.  As  a matter  of  interest,  the  generation  of  Figure  112 
required  approximately  one  hour  of  computer  time  on  a CDC  6600. 

An  interesting  next  step  for  this  analytical  effort  would  be  to 
apply  it  to  relatively  very  small  gas  turbine  engines  such  as  an  Auxiliary 
Power  Unit  (APU)  or  a Jet  Fuel  Starter  (JFS)  where  the  secondary  flows, 
leakage  flows,  etc.  are  a very  significant  portion  of  the  total  flow  rate. 
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CHAPTER  V 


CONCLUSIONS 

Several  conclusions  may  be  drawn  from  this  activity  and  the  results 
produced : 

1.  At  high  subsonic  flow  conditions,  the  Engelman  analysis  of  radial 
equilibrium  is  indeed  valid  and  such  a model  is  more  representative 
of  the  real  flow  field  than  is  the  free  vortex.  This  is  due  to  com- 
pressibility effects  and,  as  stated  by  Hawthorne  and  Horlock^'  p' 

a free  vortex  turbomachine  will  have  uniform  axial  mass  flux  distribution 
providing  the  assumption  of  incompressibility  is  valid.  Since  a free 
vortex  design  results  in  twisted  blades,  it  is  considered  significant  that 
the  nontwisted  blades  used  here  also  resulted  in  uniform  axial  mass  flux 
distribution  in  a flow  regime  where  the  assumption  of  incompressibility 
is  certainly  not  valid.  Hence,  the  conclusion  is  that  the  Engelman  pre- 
dictions were  verified  by  this  experimental  evaluation  effort. 

2.  The  Engelman  solution  is  best  characterized  as  a synthesis  tool  to 
guide  machine  designs  rather  than  an  absolute  analytical  technique.  As 
suggested  by  Hawthorne  and  Horlock^’  p'  ^0)  the  general  turbomachinery 
design  problem  "may  be  conceived  as  having  a 'direct'  and  an  'indirect' 
aspect.  In  the  former,  the  configuration  of  the  machine  is  given  and  it 
is  required  to  determine  the  behavior  of  the  fluid.  In  the  latter,  a 
particular  fluid  flow  is  specified  and  it  is  required  to  determine  the 
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machine  which  will  meet  the  specification."  Thus,  the  Engelman  tech- 
nique is  best  suited  and  intended  to  be  a basis  for  design  to  be  used 
in  much  the  same  way  as  the  basic  free  vortex  ideas  have  been.  As 
such,  it  is  applicable  to  the  so-called  'indirect'  problem. 

3.  The  overall  rig  performance  and  operational  suitability  were  quite 
acceptable.  All  features  functioned  as  designed  and  in  general  there 
were  no  significant  problems  with  taking  data  anytime  the  facility 

air  system  could  be  made  available.  The  hydraulic  servo  probe  actuator 
was  especially  convenient  and  reliable.  The  only  problem  experienced 
with  it  was  air  entrainment  in  the  fluid  after  any  long  period  of 
inactivity.  Bleeding  the  lines  always  solved  this  problem  and  caused 
no  real  hardship.  The  data  produced  by  the  test  assembly  and  its 
instrumentation  were  demonstrated  to  be  reoeatable  and  consistent. 

( 26 ) 

4.  The  two-dimensional  stack  solution  technique  developed  by  Wysong 
produces  flow  angle  predictions  whicn  are  consistently  several  degrees 
too  high  and  additionally  are  ordered  in  reverse  to  the  measured  data. 

At  these  flow  Mach  number  conditions  it  appears  that  the  combination 

of  two-dimensional  assumptions  plus  the  lack  of  allowance  for  viscosity 
effects  introduces  significant  error. 

5.  The  three-dimensional  viscous  technique  developed  by  Dodge^^ 
generates  predictions  of  flow  angle  distribution  which  are  generally 
ordered  correctly  compared  to  measured  results  but  are  still  several 
degrees  higher.  While  it  is  difficult  to  conclude  that  the  three-dimen- 
sional effort  produces  results  very  close  to  measurements  on  the  basis 
of  gross  flow  angle  comparisons,  it  does  produce  a detailed  visualiza- 
tion of  secondary  flow  effects  and  blade  loss  predictions.  In  addition. 
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the  end  purpose  of  this  three-dimensional  nroaram  is  a heat  transfer  cal- 
culation which  has  been  proven  in  other  test  cases  to  be  quite  useful. 

Due  to  the  relative  accuracy  and  completeness  of  the  three-dimensional 


analytical  program,  it  is  considered  likely  that  small  gas  turbines 
such  as  auxiliary  power  units  could  also  benefit  from  such  treatment, 
and  it  is  expected  that  such  effort  will  evolve  on  future  research 
programs  in  that  area. 
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APPEiiCIX  A 

DERIVATIONS  AMD  SOLUTIONS 

DERIVATION  OF  THE  BASIC  FOUATION 


Consider  an  annulus  utilizing  cylindrical  coordinates  as  follows: 


Then,  after  making  the  following  assumptions,  the  pressure  forces  acting 
on  differential  element  A may  be  calculated: 

(i)  Compressible,  nonviscuous  fluid. 

(ii)  Only  a combined  rotational  and  axial  motion  exist;  i.e., 
there  is  no  radial  flow. 

(iii)  The  effects  of  gravity  and  friction  are  neglected. 

(iv)  Flow  is  syrmetric  about  the  z-coordinate  axis. 

(v)  Steady  flow  throughout  the  annulus. 

( vi )  Sin  (d  ) 2:  d 
r-di rection : 

Fr  = £p  + -P-  drj  Tr  + drj  d dz  - p r d dz 

" [p  + ^r  t]  d dr  dz  (A- 1 ) 

Fr  = [r  dr  + ? |r  dr?]  d dz  (A>2) 

2 

Neglecting  the  dr  term. 

F = r dr  d dz  (A. 3) 

r c<r 

-direction: 


(A. 2) 


(A. 3) 


F ■ <*•«> 


F = o 


z-direction: 


-[«-H  dT 

fz  ' I?  iz[2  r dr  * dr?K 


(A. 5) 


(A. 6) 


(A. 7) 


By  neglecting  the  dr  term. 
Fz  = r dr  d0  dz 


(A. 8) 


The  mass  of  the  differential  element  of  fluid  is  found  from: 
dm  = density  x volume 


dm  = pdz  P it ( r + dr)2  - nr2  ] ~ 
[ 2 r dr  + dr2] 


dm  = 


pdz  do 


(A. 9) 
(A. 10) 


Neglecting  the  dr^  term. 

dm  = pr  dr  d dz  (A. 11 ) 

Acceleration  Equations 

In  cylindrical  coordinates,  the  components  of  acceleration  may  be 
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a0  = 


Then  using  the  previous  assumptions,  the  components  of  acceleration  may 
be  reduced  to:  o 


ar  = * T 


ae  = 0 


az  = 0 


(A. 15) 

(A. 16) 
(A. 17) 


Equations  of  Motion 
F = - mass  x acceleration 


r-direction: 


Fr  = - dm  ar 


(A. 18) 
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r -r^  dr  do  dz  = or  dr  dz  | - -p  | 
V2 

3£  _ p 0 

3r  g r 
ac 


C-d i recti on: 


Fe  = ' dm  ao 


o = o 


z-di rection : 


(A. 19) 


(A. 20) 


(A. 21) 


- dm  a 


(A. 22) 


r dr  d dz  = o 


(A. 23) 


Finally,  the  partial  differential  equations  of  motion  can  be  reduced  to 
a single  ordinary  differential  equation  - the  so-called  "basic"  equa- 


tion of  interest  to  this  experiment. 

, a AJL 

r gc  . 2 

pV0 


(A. 24) 


A SOLUTION  TO  THE  3ASIC  EQUATION  FOLLOWING  ENGELMAN'S  ANALYSIS 

by  making  the  following  additional  assumptions: 

(i)  The  flow  has  constant  total  enthalpy  and  constant  entropy 

throughout  the  annulus. 

(ii)  The  fluid  is  a perfect  gas 

(iii)  There  is  a continuity  of  mass 

Rearranging  the  geometric  relationship 

= + y?  (A. 25) 

0 Z 
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Where 


V2  = axial  velocity 


\L  = tangential  velocity 


V = whirl  velocity 

and  substituting  this  into  equation  (A. 24) 


dr 

- = 9C 


dp 

Dr  - M 


(A. 26) 


For  a perfect  gas  with  isentropic  flow,  the  compressible  Bernoulli 
equation  results  in: 


^ ■ 2 9c  (rr) 


Po 


1 -ffi- 

L Po. 


(¥) 


c 

°0 


(A. 27} 


Substituting  the  isentropic  relationship 

,1/k 


— - (M 

Dc  V Po  / 

into  equation  (A. 27)  and  rearranging 


1/k 


(A. 28) 


c',2-29c  (pt)po  [P,/k-p] 


(A. 29) 


Evaluating  the  geometric  relationship  from  Figure  (A-l) 

Vz  = V cos  i (A. 30) 

at  the  condition  M = 1 and  substituting  this  into  the  continuity  equation 


pVz  = p*  V 


z* 


the  following  equation  is  obtained: 


oVz  = p*  V*  cos 


(A. 31) 


(A. 32) 
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• sir  m 


(ft) 


(A. 40) 


Integrating  equation  (A. 38) from  the  condition  M = 1 to  a desired  value 

r P 

/ 

f (P)  dP  (A. 41) 

P- 


r r 

/?■/' 


Evaluating  and  rearranging 


— = exp 
■ * 


r 

I 


f (P)  dP 


(A. 42) 


The  flow  angle  can  be  expressed  by: 

-1 


$ = tan 


h) 

vz/ 


(A. 43) 


Equating  equation  (A. 24)  and  equation  (A. 38) 

9c  f (P) 

2 ' 

- pV  P0 

Evaluating  equation  (A. 43)  at  the  condition  M = 1,  substituting  into 
equation  (A.33)and  rearranging 

o 

v.  a : 

p 


(A. 44) 


pV2  . . ii  % 

° tan2** 


(A. 45) 


Rearranging  equation  (A. 43)  , equating  with  equation  (A. 45)  and 
rearranging 

2 


oV 


0 


f p*  V tan  <f.  "I 
9* 

C 

p tan  f>* 


(A. 46) 
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Dividing  gc  by  the  above  equation  and  equating  with  equation  (A. 44) 


f (P)  p tan 

Pn  = 9C 


p*  tan  * 


Substituting  the  geometric  relationship  from  Figure  (A-l) 


V = V sin  <p 


and  the  isentropic  relationship 


v = m ^ 9c  k y 

evaluated  at  the  condition  M = 1 into  equation  (A. 47)  gives 


(A. 47) 


(A. 48) 


(A. 49) 


f (P)  = p tan2  ^ 

p0  Z „ - 2 7175 


(A. 50) 


k . *.p*sin^  tan^  $ 

Substituting  equation  (A. 28),  (A. 37)  , and  (A. 39)  into  equation  (A. 50)  , 

rearranging  and  solving  for  t yields  the  following  relationship: 

1/2 

.,lc,  p’,k  I 


i = tan 


(A. 51) 


L cos  f(P)J 

Thus,  Engelman  was  able  to  determine  that  both  the  relative  radius 
ratio,  r/r#,  and  the  flow  angle,  t,  are  functions  only  of  the  pressure 
ratio,  P,  for  a given  and  k. 

A simpler  expression  relating  j#,  the  flow  angle,  i,  and  the  Mach 
number,  M,  can  be  derived  by  applying  the  previous  assumptions. 


cos  M k+1 

cos  J 2 + (k-1)  M2 


k+1 

?Tk-T) 


(A. 52) 


This  equation  was  derived  previously  in  Chapter  III. 
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Assuming  k = 1.4  and  solving  equation  (A. 52)  to  obtain: 


And  after  combining  equation  (A. 53)  and  equation  (A. 54)  to  obtain: 
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By  taking  the  following  equation  of  the  isentropic  relationship 


l/k 


= P 


l/k 


(A  .61  ) 


and  evaluating  it  at  M = 1,  the  following  relationships  can  be  obtained: 

l/k 


p* 


pT/k 


(A.62) 


dp  _ dP 
P*  " 


p; 


(A. 63) 


Then  by  substituting  the  above  two  equations  and  (A.60)into  equation  (A. 24) 
and  rearranging 


dr 


P* 


(¥) 


dP 


,2-2  l/k 
kn*  sin  : * p 


(A.  64) 


Substituting  the  isentropic  relationship 

r 2 (ft) 

£ * Ur]  ’ P* 

into  equation  (A.66)and  rearranging 


(A. 65) 


dr 

1 


k+1 


5 72 

2 kr;  sin  4>*  P 


dP 

T7k 


Integrating  from  the  condition  M = 1 to  a desired  value 

P 


/ ? = /krl  sin  7.J 

r.  P* 


dP 

„T7k 


(A. 66) 


(A. 67) 
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Integrating  the  above  equation,  substituting  equation  (A. 65),  and 
rearrangi ng . 


(¥) 


(A. 68) 


] _ (k+D  P \ 

( k-I ) sin2 


The  flow  angle  can  be  expressed  as: 


$ = sin 


’(r 


(A. 69) 


Substituting  the  isentropic  relationship 


k * [-(^)«#)  - 


(A.  70) 


along  with  (A. 57,  (A. 60),  and  (A.61)into  equation  (A. 69)  and  rearranging 


. -1 

: = sin 


» JJ3. 

r 1 P* • 

l¥hT77 


sin  f* 


(A. 71  ) 


Substituting  equation  (A. 65)  into  the  above  equation  and  rearranging 


. -1 
s i n 


sin  c* 


1 - P 


(A. 7 2 ) 
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Finally,  continuing  equation  (A. 68)  and  equation  (A. 72)  for  k - 1 .4  we 
get: 
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APPENDIX  B 


PROBE  CHECKOUT  PROCEDURE 


Even  though  a quality  wedge  probe  had  been  purchased  from  a well- 
known  manufacturer.  United  Sensor  and  Control  Corporation,  there  is 
always  the  possibility  of  a quality  control  error  which  could  affect 
the  indicated  results.  Since  the  flow  angle  was  to  be  measured  directly 
by  nulling  the  wedge  face  pressures,  it  was  considered  worthwhile  to 
set  up  some  sort  of  checkout  fixture  which  could  verify  the  accuracy  of 
such  pressure  indications.  The  plan  was  to  establish  an  axial,  uniform 
flow  field  with  a probe  mounting  scheme  located  accurately  with  respect 
to  it.  Then,  by  making  provision  to  accurately  measure  the  probe  rota- 
tion required  to  null,  a calibration  factor  could  be  established  and 
applied  to  any  indicated  flow  angle  on  the  rig. 

The  apparatus  constructed  is  shown  in  Figures  B-l,  and  B-?. 

Since  dial  indicators  were  to  be  used  against  the  face  of  the  probe 
output  head,  the  location  of  that  face  relative  to  the  wedge  faces  had 
to  be  established  accurately.  The  setup  in  the  Air  Force  Aero  Propulsion 
Laboratory  machine  shop  to  establish  these  relative  locations  using  an 
optical  dividing  head  is  shown  in  Figures  B-3,  and  B-4.  For  the  probe 
used  during  all  final  data-taking  runs,  the  output  head  face  was  deter- 
mined to  be  exactly  perpendicular  to  the  centerline  of  the  wedge  itself 
so  that  no  correction  for  that  difference  was  necessary. 
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Figure  B-l  Probe  Checkout  Fixture  - Front  View 


t Fixture  * Rear  View 


Figure  B-4  Optical  Dividing  Head  Probe  Manifold  Checkout 


The  final  operational  setup  is  shown  in  Figure  B-5.  As  can  be  seen 
there,  the  apparatus  was  connected  to  the  same  altitude  chamber  used  for 


! 


the  test  apparatus  receiver.  An  atmospheric  inlet  was  provided  to  a 
plexiglass  tube  with  an  inside  diameter  of  3 1/2  inches.  The  probe  was 
located  16  inches  or  approximately  4 1/2  tube  diameters  from  the  inlet  to 
assure  establishment  of  uniform  axial  flow.  The  probe  actuator  was 
operated  to  center  the  probe  in  the  tube  and,  of  course,  the  actuator 
mounting  pad  was  machined  very  accurately  to  assure  probe  location  relative 
to  the  flow  field.  In  addition,  the  whole  assembly, including  the  actuator, 
was  checked  by  gage  block  and  surface  plate  reference  procedures  to  assure 
the  relative  mechanical  alignment  of  the  probe  and  the  flow  field  so  that 
the  dial  indicator  readings  used  to  measure  rotation  would  be  direct 
answers . 

The  end  result  was  that  the  probe  used  for  all  runs  in  this  disserta- 
tion was  perfect  within  the  capability  to  measure  error  described  above. 
Therefore,  no  correction  was  applied  to  the  indicated  results.  It  is 
worth  noting,  in  closing,  that  the  spare  probe  had  an  error  of  approximately 
3°  using  this  same  checkout  procedure. 


i 
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Figure  B-5  Probe  Checkout  Flow  Setup 


APPENDIX  C 


ACTUATOR  CALIBRATION 

To  make  sure  the  probe  would  come  as  close  to  the  centerbody  as 
possible  without  damaging  contact,  a piece  of  0.020"  shim  stock  was 
taped  to  the  centerbody  as  a spacer.  The  probe  was  then  moved  circum- 
ferentially to  determine  the  highest  point  and  the  probe  was  then 
scribed  to  make  sure  it  could  not  be  set  lower  than  the  highest  point 
encountered.  This  reference  mark  was  then  used  to  set  the  probe  at  the 
correct  height  for  all  calibrations  and  all  installations  in  the  rig 
itself. 

The  angular  calibration  of  the  probe  was  conducted  by  placing  a 
round  calibration  disk  (degree  wheel)  on  the  probe  itself.  This  can  be 
seen  in  Figure  C-l.  This  disk  was  graduated  in  degrees  and  by  use  of  the 
vernier,  readings  could  be  taken  in  0.1  degree  increments.  The  probe 
was  hand  turned  relative  to  the  actuator  to  read  50  degrees  relative  to 
the  centerbody  center  line  which  would  now  be  the  lowest  degree  reading 
possible;  i.  e.  the  starting  point  of  the  rotary  actuator  travel  range. 
Calibration  was  then  started.  One  operator  was  stationed  looking 
through  a magnifying  lamp  at  the  disk  and  another  turned  the  probe  con- 
trol until  told  to  stop  by  the  first  operator.  If  the  movement  was  too 
small  or  too  large,  iterative  adjustments  were  made  until  the  desired 
setting  was  obtained.  The  digital  voltmeter  reading  was  then  recorded. 
This  continued  for  every  calibration  point  required.  The  next  day,  the 
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second  operator  positioned  the  probe  by  use  of  the  voltmeter  and  the 
first  operator  observed  the  degree  wheel  reading,  therefore  checking  the 
calibration.  No  significant  differences  were  encountered  and  the  read- 

N 

ings  were  accepted.  The  calibration  readings  taken  here  and  used  for 
all  tests  are  shown  in  Table  C-l . 

The  linear  calibration  was  done  in  the  manner  shown  in  Figure  C-2 
One  operator  looked  through  a magnifying  lamp  and  told  a second  operator 
when  the  dial  moved  1/10  of  an  inch.  The  voltmeter  reading  was  then 
recorded.  This  was  done  for  22  readings  and  a final  reading  which  was 
not  1/10  of  an  inch.  These  readings  were  all  checked  in  the  same  manner 
used  to  check  the  angular  calibration.  They  are  also  shown  in  Table  C-l. 

p. 


Table  C-l 

ACTUATOR  CALIBRATION 


LINEAR 

ANGULAR 

VOLTMETER 

(VOLTS) 

POSITION 

(INCHES) 

VOLTMETER 

(VOLTS) 

POSITION 

(DEGREES) 

1.956 

0.29 

0.315 

50 

2.193 

0.39 

0.400 

51 

2.444 

0.49 

0.493 

52 

2.701 

0.59 

0.574 

53 

2.966 

0.69 

0.687 

54 

3.218 

0.79 

0.786 

55 

3.448 

0.89 

0.859 

56 

3.886 

0.99 

0.958 

57 

3.916 

1.09 

1.069 

58 

4.160 

1.19 

1.150 

59 

4.398 

1.29 

1.248 

60 

4.642 

1.39 

1.354 

61 

4.P94 

1.49 

1.415 

62 

5.133 

1.59 

1.536 

63 

5.378 

1.69 

1.580 

64 

5.625 

1.79 

1.691 

65 

5.866 

1.89 

1.760 

66 

6.113 

1.99 

1.860 

67 

6.358 

2.09 

1.953 

68 

6.627 

2.19 

2.024 

69 

6.857 

2.29 

2.129 

70 

7.066 

2.39 

2.226 

71 

7.227 

2.45 

2 280 

72 
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APPENDIX  D 

NUMERICAL  INDEX  OF  TEST  APPARATUS  DRAWINGS 


r 

! ' 

t 


DWG.  No 

Date 

Title 

7552-04-2663 

5 

Jun 

74 

Blade  Ass'y  Installation 

2664 

5 

Jun 

74 

Plate-Hold  down 

2665 

5 

Jun 

74 

Spacer 

2666 

10 

Jun 

74 

Probe  Drive  Ass'y 

2667 

17 

May 

74 

Drive  Saddle 

2668 

24 

May 

74 

Bar  Screw  Ass'y 

2669 

5 

Jun 

74 

Plate-Actuator  Mounting 

2670 

3 

Jun 

74 

Motor  Mount 

2671 

21 

May 

74 

Angle-Roller  Support 

2672 

24 

May 

74 

Block-Roller  Support 

2673 

3 

Jun 

74 

Drive  Pin 

2674 

3 

Jun 

74 

Ball  Bushing  Pushrod 

2675 

28 

May 

74 

Bracket 

2676 

7 

Jun 

74 

Pointer 

2677 

7 

Jun 

74 

Spacer-Pointer 

2678 

29 

May 

74 

Ball  Bushing  Housing 

2679 

30 

May 

74 

Ball  Bushing  Shaft  Support 

2680 

21 

May 

74 

Plate  Ass'y  Cylinder  Lock 

2681 

10 

Jun 

74 

Degree  Rule 

2682 

5 

Jun 

74 

Plug-Clamp 

2683 

5 

Jun 

74 

Clamp 

2684 

5 

Jun 

74 

Stud 

2685 

5 

Jun 

74 

Nose  Cone 

2686 

5 

Jun 

74 

Blade  Ass'y 

2687 

5 

Jun 

74 

Blade  Ass'y  69.7° 

2688 

5 

Jun 

74 

Blade  Ass'y  65° 

2689 

5 

Jun 

74 

Blade  Ass'y  60° 

2690 

5 

Jun 

74 

Ring,  Bellrnouth  Mounting 

2691 

5 

Jun 

74 

Ring 

2692 

5 

Jun 

74 

Liner 

2693 

5 

Jun 

74 

Shell 

2694 

5 

Jun 

74 

Cyl inder 

2695 

5 

Jun 

74 

Tube-Support 

2696 

5 

Jun 

74 

Pressure  Port  Ass 'y 

2697 

5 

Jun 

74 

Support  Ass 'y  Aft. 

2698 

5 

Jun 

74 

Mounting  Plate 

2699 

5 

Jun 

74 

Mounting  Saddle 

2860 

3 

Sep 

74 

Motor  Support  Plate 

2861 

3 

Sep 

74 

Motor  Mounting  Plate 

2862 

3 

Sep 

74 

Motor  Support  Angle 

2863 

4 

Sep 

74 

Ass'y  Bar  Screw 

2864 

4 

Sep 

74 

Support  Angle 

- . — LdL. 
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APPENDIX  E 


Hi  101+  o Cj«  157?.  tt  &■>  n.V  i$(.  mv  irii  rfo  W*  rJ-.5>  '/.tf  If  /*? 

17.7.  it.  f .0.  Kp  isrt.  n ins  irt.i  ifi,j  jmtn, 1 it>>  iym  iv  i}.«  /U 

ft"/  V->  .0.  057  |.$>7  70  ff  iq.1  is$.f  wil  If%i  ifp  /It-I  I'fl-'l  . l)fi  ll-K  °lu  III 

^■7,  f.'  l .(?.  2chJ4  '.i !£-.  7°  ??  ij»°  'rtf  iff*  iris  isis  mi  /y?A  . >5°  1 s«,rfV 

IIES  1i-t  ioi  7VJ  i-?7^  70  rt  1*1  isi./  1 fit  If,,y  [fri  (fa  up  111  ,1>  ly  rv% 


Table  1-2  0°  - Probe  on  Vertical 


APPENDIX  F 


LONG  VS  SHORT  TRAVERSE  COMPARATIVE  PLOTS 

Presented  in  this  appendix  are  plots  of  data  taken  at  station  three 
using  the  69.7°  cascade.  The  circumferential  location  was  2.5  degrees 
right  of  the  reference  blade  trailing  edge  and,  for  comparison,  both 
"long"  and  "short"  traverses  were  made.  The  radial  positions  included 
in  these  traverses  was  shown  previously  as  Table  4 in  Chapter  II. 

As  can  be  seen  in  Figures  F-l  through  F-3,  the  short  traverse  pro- 
duces data  which  is  quite  comparable  to  that  produced  by  the  long 
traverse.  Hence,  in  order  to  conserve  energy  and  generate  the  maximum 
results  in  the  minimum  time,  much  of  the  data  produced  during  this 
experiment  was  taken  with  a short  traverse. 


Blade  Angle  69.7° 
Station  3 


r 


PROBE  POSITION  (IN  PROM  CENTERBODY) 


Blade  Angle  69.7° 
Station  3 
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APPENDIX  G 


TOTAL  PRESSURE  AND  FLOW  ANGLE  PLOTS 
FOR  CHECKOUT  RUNS 

The  data  shown  here  was  taken  during  shakedown  runs  on  the  test  appar- 
atus and  was  used  to  optimize  the  size  of  backflow  plates  (flow  fences) 
for  the  channel  exit  plane  for  each  cascade.  These  fences  were  attached 
to  the  wood  liner  at  the  test  apparatus  exit  plane  and  cut  down  the 
channel  outside  diameter  to  the  particular  size  of  the  backflow  plate  bore. 
Figures  G-l  and  G-2  show  total  pressure  and  flow  angle  respectively 
across  the  channel  for  varying  exit  diameters.  These  data  are  all  for 
the  69.7°  cascade.  Shown  in  Figure  G-3  is  an  enlarged  plot  of  the 
selected  plate  bore  (14.625  inches)  results;  i.e.  total  pressure  varia- 
tion across  the  channel.  In  this  case,  for  the  intended  checkout  purpose, 
only  raw  data  was  plotted  so  that  the  ordinate  in  volts  was  taken 
directly  from  the  probe  actuator  readout  without  conversion  to  actual 
geometric  channel  position. 

Finally,  in  Figure  G-4  is  shown  the  comparison  of  effects  of  no 
backflow  plate  vs  the  selected  one  at  15.403  inches  on  total  pressure 
profile  for  the  65°  cascade. 

Further  checkout  runs  confirmed  the  trend  shown  by  these  results; 
i.e.  no  plate  was  required  to  produce  a very  constant  total  pressure 
profile  across  the  channel  when  the  6C°  cascade  was  installed. 
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Figure  G-l  Total  Pressure  VS  Probe  Position  - 69.7°  Cascade 


Shown  in  Figure  G-5  is  a typical  flow  fence  made  of  wood, while  in 


Figure  G-6  one  made  of  aluminum  is  shown  installed  in  the  rig.  The 
"14.25"  stamped  into  the  wood  is  the  inside  diameter  of  the  fence.  The 
original  plan  was  to  use  wooden  ones  while  sizing  the  exit  area  and  then 
make  aluminum  ones  which  would  be  durable  enough  for  unlimited  running. 

However,  experience  proved  that  the  wooden  ones  were  quite  durable  and 
they  were  used  for  most  runs. 

i 

t 


Figure  G-6  Flow  Fence  Installation 
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appendix  h 

SURFACE  ROUGHNESS  MEASUREMENTS 


Measurements  of  surface  roughness  were  made  on  representative 
sections  of  all  parts  of  the  rig  flowpath.  The  measurement  locations 
and  results  are  shown  as  follows: 


Liner 


Centerbody 


Figure  H-l 
Figure  H-2 


69.7°  Blades 


Table  H-l 


65°  Blades 


Table  H-2 


60°  Blades 


Table  H-3 


• t / • 


Table  H-2 


Surface  Roughness  - 65#  Blades 


BLADE  NUMBER  1 


DATA  SETS 

READING  POSITION  (measured  in 

RMS)* 

1 

2 

3 

4 

1 

45 

70 

■K9| 

— 

2 

40 

65 

3 

40 

65 

HH 

IfH 

BLADE  NUM3ER  17 


DATA  SETS 

READING  POSITION  (measured  in  RMS)# 

1 

2 

3 

1 

65 

80 

2 

70 

90 

Ifa 

3 

65 

135 

80 

Ha 

BLADE  NUMBER  33 


DATA  SETS 

READING  POSITION  (measured  in 

RMS)# 

1 

2 

3 

4 

1 

H 

90 

80 

2 

til 

95 

100 

3 

la 

100 

90 

120 

All  surface  roughness  measurements  are  in  microinches. 
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APPENDIX  I 


REYNOLD'S  NUMBER  CALCULATIONS 

Reynold's  numbers  based  on  twice  the  channel  height  (twice  the 
blade  span)  were  calculated  over  the  entire  range  of  flow  rates  inves- 
tigated. Table  1-1  summarizes  these  results  and  the  actual  calculations 
follow.  In  the  first  set  of  calculations  the  extreme  cases  of  velocities 
achieved,  both  high  and  low,  were  used  and  in  the  second  set  average 
values  were  used.  The  velocities  were,  of  course,  determined  from  the 
wedge  probe  pressure  data  and  are  in  the  measured  direction  of  whirling 
flow.  The  sunnary  table  is  based  on  the  second  set  of  calculations  and 
the  first  set  is  included  here  for  extreme  bounding  information  only. 


Table  1-1 

Flow  Ranges  Investigated 


Reynold's  Number  Calculations  - Extreme  Cases 


60°  Blades  Low  Flow 

Re  = (698  ft/sec)  (0.071  lbm/ft3)  (5/12  ft)  = -|  607  x 106 

1.285  x 10"^  lbm/ft  sec 

Re  = (383^)_lfh0711J5Z12l  s g g3  x 1Q5 

1.285  x 10‘5 


65°  Blades  Low  Flow 

Re  = (416.5)  (0.071)  (5/12)  = g5g  x }QS 

1.285  x 10'5 

Re  = (648.45)  (0.071)  (5/12)  = ! 493  x 1Q6 

1.285  x 10‘5 


65°  Blades  Low  Flow 

Re  = (473.24)  (0,071)  (5/12)  = ] Q8g  x ]Q6 

1.285  x 10'5 

Re  - [719.9)  (0.07TJL(5/X2i  . ,.657  x 106 

1.285  x 10'b 


60°  Blades  High  Flow 

Re  - (494,2)  (0.071)  (5/12)  . 1J38  x 1Q6 

1.285  x 10"b 

Re  = (701.9)  (0.071)  (5/12)  = t 616  x 1Q6 

1.285  x 10"5 


65°  Blades  High  Flow 


Re  = (465)  (0.071)  (5/12)  = 1-07  x 106 

1 .285  x 10'5 

Re  = (844.51)  (0.071)  (5/12)  = 1944  x 1Q6 

1.285  x 10'5 


69.7°  Blades  High  Flow 
Re  = (701.26)  (0.071)  (5/121  = 1614  x 1Q6 


-5 


1.285  x 10 

^ = (828.9)  (0.071)  (5/12)  = 1 .908  X 1 06 

1.285  x 1C'5 

Reynold's  Number  Calculations  - Average  Cases 
60°  Blades 


Re  = (680  ft/sec)  (.071  lbm/ftJ)  (5 

1.285  x 10"5  lbm/ft-sec 


= 1.5655  x 10 


Re  = (500  ft/sec)  (.071  lbm/^t3)  (5/12  ft)  = 1 15]1  x 1Q6 

1.285  x 10*5  lbm,  ft-sec 

1.5655  x 106 
Range  - Re  = ] .1511  x IQ6 


65°  Blades 

Re  = (830)  (.071  lbm/ft3)  (5/12  ft)  . ] 91Q8  x 1()6 

1.285  x 10-5  lbm/ft-sec 

Re  = (570  ft/sec)  (.071  lbm/ft3)  (5/12  ft)  = ] 3123  x 1()6 

1.285  x 10'5  lbn/ft-sec 

1 .9108  x 106 
Range  - Re  * c 

1 .3123  x 10° 
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i 69.7°  Blades 

Re  = (850  ft/secU_.Q71  lbm/ft3)  (5/12  ft)  = i.9569  x -|06 
1.285  x 10'5  Ibm/ft-sec 

Re  = (615  ft/sec)  (.071  lbm/ft3)  (5/12  ft)  = •,  41586  x 106 
1.285  x 10‘5  Ibm/ft-sec 

1 .9569  x 106 

Range  -*•  Re  = 

1.41586  x 10° 


APPENDIX  J 


MASS  FLOW  CALCULATION  PROCEDURE 


A computer  program  was  written  which  calculated  the  mass  flow  as  a 

function  of  probe  position  and  then  summed  to  get  the  total  mass  flow  in 

the  channel.  While  calculating  the  mass  flow,  other  useful  information 

was  generated  and  printed  for  each  case  as  follows:  Measured  Mach  num- 

3 

ber.  Axial  Mach  number  and  velocity  (ft/sec).  Density  (lb/ft  ),  and  Pres- 
sure ratio  PS/PT. 

The  information  input  to  the  computer  was  the  area  slice,  total 
pressure,  pressure  difference,  flow  temperature  and  flow  angle.  The  com- 
plete program  is  listed  in  Table  J-l , and  the  line  numbers  referred  to 
in  the  following  description  of  operation  are  shown  there.  The  program 
first  forms  the  pressure  ratio,  (lines  14,  27,  and  28), 


Pressure  Ratio  = PCC 


and  from  this  calculates  the  Mach  number  (line  29), 


Mach  Number  = A 


■1/y  1/2 

-1  x5 


.2857 


Since  for  mass  flow  purposes,  the  axial  velocity  is  desired,  first 
the  axial  Mach  number  is  found  (line  30), 

Axial  Mach  Number  = AX  = A cos  C 

and  then  the  axial  velocity  is  computed  (lines  25  and  31), 


• • • 


• % * 


• • 


o • 


Axial  Velocity  = A cos  e 


JyRgcTE 


= AxJl  .4(53.4)  (32.2)  TE 

The  density  is  then  calculated  (line  33)  assuming  an  ideal  gas, 

P = oRT 


or 


D 


P PS (5 . 1 98) 

RT  - 53T4TTE) 


= RHO 


Mass  flow  rate  (line  35)  can  now  be  calculated, 

mass  flow  rate  = FLOW  = cAV  = RH0(AR)V 
Flow  rate/area  is  calculated  by  line  37  , 

Flow  rate/area  = pAV/A  = FLOW/AR  = ARF 
The  total  mass  flow  then  comes  from  adding  the  various  contributions 
from  the  individual  slices  of  area  and  is  done  on  line  39, 

Total  mass  flow  = FLOWS  - (FLOW),  + (FL0W)2  + . . . + (FL0W)?3 


The  remaining  cards  transform  the  various  parameters  into  proper  units  or 
introduce  various  constants  or  are  logic  and  control  cards.  The  program 
can  handle  N sets  of  date  that  consist  of  23  data  points. 
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Table  J-l 


Flow  Rate  Calculation  Program 


5 


10 


15 


20 


25 


30 


35 


40 


50 


55 


61 


65 


PROGrmM  It,  IP (INPUT  .OUT3 JT , T A»t 5 = 1 NPUT ,T A PE  6 = OUT PUT ) 

k£AD(  5, 12)  N 

J=1 

51  F.EA0(5,6)SA,SB,SC,SC,SE,  SF,SG 
WRII  t(6,  HISA,S3,SC,SO,St 
WrIT  F (6,19)  SF  ,SG 
FLOWS=0 

C WRITE  T He  HL  AGINGS 

NKITl  ( 6,21) 

WRITE (6,23) 

1=1 

I REAP ( 5,lu) Ah, FT,Pu, T ,TH£ TA 

C READ  THE  ARtA,  TEMP.,  AsO  THE  FLOW  ANGLE  FOk  THE  CALCULATIONS 

PS  = FT  -FO 

C CONSTANTS 

GAM  = 1. 4 
GA  S=  53, 4 
GKA=  32,2 

C CONVERSION  OF  1 TO  DEGREES  R 

TE=T ♦459,67 

C CONVERSION  FACTOR 

7TMETA=THETA*3.141592694/150. 

C TTHET*  IS  RADS 

OTHLTA=COS(TTHETA) 

ROOT=SQRT tGAM*GAS»ikA»TE ) 

C CALCULATION  OF  VELO- ITY - AXIAL 

PC= ( PT /PS ) 

PCO= 1/PC 

A = SORT ( ( <PC**.2857>-1>,5) 

A X = A*CT  HE  T A 
V=A*R00T»CTHETA 

C CALCULATION  OF  jENsITY 

RHO=PS»  5.  190/  (GAS»tl  ) 

C CALCULATION  OF  NASs  FLOW  RATE 

flow=rho*lr*v 

C CALCULATION  OF  FLOW  RATE/  UNIT  AREA 

A„F=F  LOW/AR 

C TOTAL  MASS  FLOW  RA1ES 

FLOWS=FLOWS*FLOW 

WRITE (6, 22' AR.THET  A, A ,AX ,RHD, FLOW, ARF , V,PCC 
IFCI.tQ.23)  GO  TO  999 
I = I»1 
GO  TO  1 

999  WRITE  (6,24)FLCWS 

200  IF  ( J.  FQ.  N)  GOT  0 20 1 
J=  J«  1 
GO  TO  51 
2Q1  STOP 

10  FORMAT  ( 5f  10.6) 

12  FORMAT  (12) 

15  FORMAT  <"  X NOT  CAlIBRATEO  FOR  THIS  NUMBER  ",F10.5> 

8 FORMAT  (312,F5.2,I2,F5.2,F4.1) 

II  FORMAT  ( 1H, ///, 9X, "DATf, 2X, 12,"/", 12,"/“, 12, 5X,"  POSITION*  ",F5.2, 
1"  OE&REtS  LEFT  OF  ZERO", 5X,"SAUuLE  STATION  ",I2) 

19  FORMAT (28x, "CHAH8ER  PRESSURE  ",F5.2,"  IN  HG" , 5X, "8L AO E ANGLE  ", 
1F5.2,"  DEGREES") 

21  FORMAT ( 1H,///,9X, -Area*,  4X, "FLOW  ANGL t", 2X , "MACH  NUMBER-", 2X, 
l-MACM  NUMBER-*, 2 X , " JENil T Y", IX, "FLOW  RAT E" , 3X , "FL ON  R AT  E/AREA" , 2X , 
1" VELOCITY- A XI AL",2X, “PS/ PT") 

24  FORMAT  (9X, "TOTAL  FLOW  RATt  IS  ",F12.6,"  LB/SEC") 

22  FORMAT  C8X,F7.6,2X,F5.2,7X,F10.8,  4X,F10.8,4X,F7.6,3X,F8.6,4X,F11.6 
1,4X,F9.4,7X,F7.6) 

23  FORMAT C9X,"FT**2",3X, "DEGREES", 5X,"bLA0E  EX  IT", 4X , "AXIAL",  SX, 
l"Ld/FT«»3",2X,"L«/StC",6X,"L3/SEC/FT**2",4X,"FT/SEC"> 

ENO 
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APPENDIX  K 

CHECK  RUN  AT  STATION  SEVEN 

Presented  in  this  appendix  are  the  results  from  a traverse  made 
at  station  seven  for  the  purpose  of  pursuing  further  the  general  observa- 
tion that  relatively  little  flow  straightening  was  occurring  between 
stations  one  and  three.  As  discussed  in  Chapter  IV,  these  results  do 
indeed  verify  that  relatively  little  flow  straightening  occurs  in 
several  chord  lengths  downstream  from  the  cascade. 

l 
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Blade  Angle  65.0° 
Station  7 


Figure  K-l  Ptota1  VS  Probe  Position 
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PROBFZ  POSITION  (IN  PROM  CENTERBODY) 


Blade  Angle  65.0° 
Station  7 
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Blade  Angle  65.0° 
Station  7 


50  0 52  0 54-0  56-0  58  0 BOO 


FLOW  ANGLE  (DEG) 

Figure  K-3  Flow  Angle  VS  Probe  Position 


220 


lable  K-l 


APPENDIX  L 


ERROR  ANALYSIS 

All  pressure,  temperature  and  flow  angle  measurements  in  this  experi- 
ment were  set  up  in  accordance  with  the  data  and  recommendations  reported 
f?4) 

by  Dean  and  applied  as  closely  as  possible  to  this  particular  experi- 
mental apparatus.  Other  references  used  were  John^"®’  PP*  330-338)  ancj 

(29) 

United  Sensor  and  Control  Corporation  . A standard,  commerci ally- 

avail  able  wedge  probe  manufactured  by  United  Sensor  and  Control  Corporation 

was  chosen  as  the  traversing  pressure  and  temperature  sensor.  It  was 

identified  by  the  manufacturer  as  part  number  WT  250-24-CD-C/C-36  and  was 

configured  with  a wedge  at  the  end  of  a tube  0.25  inches  in  outer  diameter. 

(29) 

Based  on  United  Sensor  and  Control  Corporation'  ' experience  and  recom- 
mendations, the  effect  of  probe  size  as  a function  of  channel  area  can  be 
determined  as  follows: 

"For  round  probes  perpendicular  to  the  flow,  the  drop  in  static 
pressure  at  the  probe  cross-section  and  downstream  from  it  is 
approximately:" 


Where:  P and  P are  the  original  and  modified  static 
S1  S2 
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pressures  respectively 


Pt  = Total  pressure 
a = Cross  sectional  area  of  probe 
A = Channel  area 

. o 

For  a probe  diameter  of  0.25  inches,  a = .04906  in  and  the  annulus  area 

A = 105.98  in^  therefore, 

J = .00046 
and 


.00056 


Hence,  the  conclusion  that  the  disturbance  effect  of  the  probe  on  static 
pressure  readings  may  be  neglected  is  reasonable. 

Based  on  the  work  of  Keast^^,  if  the  location  of  the  direction- 
measuring taps  in  the  wedge  faces  is  such  that: 


7gcV 


ps 

2 


= 0 


Where: 

Pi  = Indicated  Static  Pressure 
P$  = Actual  Static  Pressure 
p = Density 
V = Velocity 

gc  = Gravitational  Constant 

then  the  static  pressure  can  be  read  directly.  For  this  relationship  to 
be  satisfied,  Keast's  data  shows  that  the  ratio  of  the  distance  of  the  taps 
from  the  leading  edge  of  the  wedge  to  the  overall  length  of  the  wedge  face 
must  equal  0.375.  ^e^'urement  of  the  probe  used  in  this  experiment 
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indicated  that  this  relationship  was  satisfied. 

It  is  recognized  (see  John^28,  p>  331  that  the  most  ideal  static 
pressure  probe  would  have  a considerably  longer  and  conical  nose  section 
with  the  measuring  tap  located  10  to  20  probe  diameters  downstream  from 
the  nose.  However,  such  a probe  would  be  totally  impractical  for  use  in 
this  experiment  and,  the  use  of  a wedge-type  probe  was  considered  a very 
reasonable  compromise.  The  probe  manufacturer.  United  Sensor  and  Control 
Corporation,  has  stated  in  their  catalog  that  this  probe,  TYPE  WT,  is 
usable  up  to  Mach  0.7  and  in  general,  that  value  was  not  exceeded  by  more 
than  5%  in  any  of  the  traverses  made  in  this  experiment.  However,  as  a 
final  observation,  Keast^8,  p‘  887^  reports  that  his  wedge  probe  "measures 
the  average  static  pressure  to  within  1.5  per  cent  when  used  for  a wake 
traverse  of  a turbine  cascade  at  an  outlet  Mach  number  of  0.8  and  was 

even  closer  at  lower  Mach  numbers." 

( 29 1 

The  calibration  data  bulletin'  ' supplied  by  the  manufacturer  of  the 
wedge  probe  discusses  several  sources  of  error  in  static  pressure  measure- 
ment. These  are: 

Mach  number  effects 
Pitch  Angle 
Yaw  Angle 
Immersion  Depth 

The  yaw  angle  error  is  declared  by  the  manufacturer  to  be  zero  when  the 
probe  has  been  rotated  such  that  the  wedge  face  taps  indicate  equal 
pressure.  Specific  calibration  data  was  supplied  on  Mach  number  effects 
and  pitch  angle  effects.  For  the  probe  used  in  this  experiment,  the  pitch 
angle  effects  cancelled  the  Mach  number  effects  almost  completely  in  the 
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range  of  Mach  numbers  of  greatest  interest;  i.e.  0.6  and  above.  The  remain- 
ing effect  of  immersion  depth  is  presented  as  "due  to  effects  of  boundaries 
of  the  passage  and  the  secondary  flow  along  the  upstream  edge  of  the  probe 
which  in  turr  is  influenced  by  the  total  pressure  gradient  in  the  passage." 
For  probes  cantilevered  out  from  one  side  of  a passage  (the  case  for  this 
experiment),  it  is  stated  that  the  probe  will  read  true  static  pressure 
"when  the  holes  (taps)  are  ten  probe  diameters  away  from  the  boundary 
they  project  through  and  the  probe  tip  is  at  least  one  diameter  from  the 
opposite  wall."  The  specific  data  supplied  on  that  source  of  error  shows, 
for  the  size  channel  used  in  this  experiment,  a maximum  error  of  1.5%  of 
indicated  static  pressure  for  the  range  of  static  pressures  measured. 

The  readout  of  probe  pressure  indications  was  accomplished  by  visual 
observation  of  Kollsman  aneroid  gages  graduated  in  inches  of  water.  These 
were  calibrated  by  the  Precision  Measuring  Equipment  Laboratory  in  accor- 
dance with  the  provisions  of  Technical  Order  (TO)  00-20-14,  Air  Force 
Metrology  and  Calibration  Program.  Specifically,  the  Kollsman  gages  are 
category  3 items.  The  standards  used  are  directly  traceable  to  the 
National  Bureau  of  Standards  (NBS)  through  the  Directorate  of  Metrology, 
Newark  Air  Force  Station,  Air  Force  Measurement  Standards  Laboratory,  and 
the  Wright-Patterson  Air  Force  Base  Precision  Measurement  Equipment 
Laboratory  (PMEL). 

The  inner  wall  (centerbody)  static  pressures  were  read  from  single 
tube  vertical  mercury  manometers  graduated  in  tenths  of  an  inch, and  the 
scales  were  "zeroed"  prior  to  each  run. 

The  outer  wall  static  pressures  were  read  from  Kollsman  aneroid  gages 
identical  to  those  used  for  the  wedge  probe.  Therefore,  the  accuracy 
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discussion  presented  for  the  wedge  probe  gages  is  directly  applicable  and 
will  not  be  repeated  here. 


r 


j 


Both  stream  total  temperature  and  air  inlet  temperature  were  measured 

with  thermocouples.  The  total  temperature  thermocouple  was  built  into  the 

wedge  probe  while  the  air  inlet  thermocouple  was  suspended  approximately 

5 feet  in  front  of  and  4 feet  above  the  bellmouth  centerline. 

(24  d 521 

These  were  copper-constantan  thermocouples  and  Dean  ’ K' 
reports  a standard  limit  of  error  range  of  + 1 1/2°F  for  the  temperature 
range  of  -75  to  +200°F.  Dean^’  p'  ^ classifies  all  common  temperature 
measuring  situations  into  three  ranges  as  follows: 

Range  1 : Measurements  in  low  velocity  streams  with  negligible 
radiation. 

Range  2:  Measurement  in  high  velocity  streams  with  negligible 
radiation . 

Range  3:  Measurements  under  conditions  specified  in  1 or  2 but 

with  significant  radiation  effects, 

if 

For  the  two  temperatures  required  by  this  experiment,  Range  1 covers  the 
air  inlet  measurement,  and  Range  2 covers  the  stagnation  temperature 
measurement  at  each  traverse  position. 

Measurements  within  Range  1 are  reported  to  be  achievable  with  bare 
thermocouples  or  mercury  thermometers  with  negligible  error.  In  the 
absence  of  radiation  effects,  measurements  within  Range  2 are  reported  to 
be  achievable  over  a wide  range  of  subsonic  and  supersonic  Mach  numbers 
with  a correction  factor  as  high  as  0.99.  This  is  accomplished  by  various 
designs  of  vented,  diffusion-type  containers  for  the  thermocouple  inside 
the  probe, and  the  United  Sensor  Type  WT  probe  used  in  this  experiment  fits 
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these  specifications  suitably  for  the  Mach  number  range  explored  here. 

Hence,  the  conclusion  that  correction  of  the  raw  temperature  data  is 
not  required  is  considered  reasonable. 

Calibration  of  the  entire  temperature  instrumentation  subsystem  was 
conducted  by  tne  Air  Force  Aero  Propulsion  Laboratory  Instrumentation  Shop. 
The  Brown  recorder  is  a category  2 item  and  was  calibrated  using  a PMEL- 
calibrated  32°F  reference  junction  and  a PMEL-cal ibrated  precision  potentio- 
meter. International  Practical  Temperature  Scale  of  1968  (IPTS-68)  conver- 
sion tables  were  used  to  convert  the  thermocouple  readings  to  actual 
temperature  values.  Thus,  system  traceability  is  established  back  to  NBS 
just  as  for  the  pressure  gaqes. 

Measurement  of  both  the  radial  location  of  the  probe  in  the  flow 
channel  and  the  flow  angle  by  use  of  a hydraulic  servoactuator  system  was 
selected  due  to  the  convenience  and  proven  repeatability  of  such  devices. 

The  calibration  of  this  system  in  both  linear  and  angular  modes  is 
covered  in  detail  in  Appendix  C.  In  addition,  determination  of  the 
so-called  "zero  reference"  for  the  particular  wedge  probe  used  is  pre- 
sented in  Appendix  B.  Selection  of  the  "null  method"  rather  than  the 
"fixed  method"  as  described  by  Dean^’  was  made  to  eliminate  the 

requirement  *or  any  compressibility  correction  on  flow  angle  indications 
from  the  wedge. 

In  summary,  it  is  estimated  that  the  accumulated  total  error 
on  each  parameter  is  bounded  as  follows: 

Flow  Angle  +0.25° 

Stream  Total  Pressure  1*  of  indicated  pressure 

Stream  Static  Pressure  2%  of  indicated  pressure 
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Outer  Wall  Static  Pressure  +0.5  inches  of  water 

Inner  Wall  Static  Pressure  +0.1  inch  of  mercury 

Total  Temperature  +1  1/2°F 

Air  Inlet  Temperature  +1  1/2°F 

Circumferential  Position  +1/4° 

Radial  Position  +.005  inches 

The  equations  and  computer  prediction  models  governing  the  results  and 
conclusions  presented  in  this  Dissertation  are  such  that  the  most  critical 
measurement  is  clearly  the  flow  angle.  Accurate  determination  of  flow 
angle  was  a primary  objective  from  the  outset  and, it  is  considered  reason- 
able to  conclude  that  very  acceptable  accuracy  was  achieved.  The  next 
critical  measurement  was  pressure  and  within  the  limits  of  resources 
available  for  this  experiment,  the  error  is  considered  acceptable  to 
support  the  conclusions  drawn.  The  error  bounds  established  for  the 
temperature  measurements  will  result  in  no  significant  change  to  the 
conclusions  drawn  here. 

Finally,  with  respect  to  the  circumferential  location  in  the  channel 
and  the  radial  traverse  points,  the  errors  here  have  the  least  possible 
effect  on  the  conclusions.  However,  as  presented  in  the  main  body  of  the 
Dissertation,  repeatability  of  results  was  verified  several  times  during 
the  course  of  the  experimental  phase  and  was  proven  to  be  quite  satisfac- 
tory. This  verification  required  returning  to  a given  circumferential 
location  and  repeating  a radial  traverse. 


229 


BIBLIOGRAPHY 


1.  Wu,  Chung-Hua,  Wolfenstein,  L.,  "Application  of  Radial -Equi I ibrium 
Condition  to  Axial-Flow  Compressor  and  Turbine  Design,"  NACA 
Report  955,  1950. 

2.  Wu,  Chung-Hua,  "A  General  Theory  of  Three-Dimensional  Flow  in  Sub- 
sonic and  Supersonic  Turbomachines  of  Axial-,  Radial-,  and  Mixed-Flow 
Types,"  NACA  Technical  Note  2604,  1952. 

3.  Wu,  Chung-Hua,  "A  General  Through-Flow  Theory  of  Fluid  Flow  with 
Subsonic  or  Supersonic  Velocity  in  Turbomachines  of  Arbitrary  Hub 
and  Casing  Shapes,"  NACA  Technical  Note  2302,  1951. 

4.  Stanitz,  J.  D.,  "Approximate  Design  Method  for  High-Solidity  Blade 
Elements  in  Compressors  and  Turbines,"  NACA  Technical  Note  2408,  1951. 

5.  Wu,  Chung-Hua,  Brown,  C.  A.,  "Method  of  Analysis  for  Compressible 
Flow  Past  Arbitrary  Turbomachine  Blades  on  General  Surface  of 
'Revolution,"  NACA  Technical  Note  2407,  1951. 

6.  Huppert,  M.  C.,  MacGregor,  C.,  "Comparison  Between  Predicted  and 
Observed  Performance  of  Gas-Turbine  Stator  Blade  Designed  for 
Free-Vortex  Flow,"  NACA  Technical  Note  1810,  1949. 

7.  Smith,  L.  H.,  Jr.,  Traugott,  S.  C.,  Wislicenus,  G.  G.,  "A  Practical 
Solution  of  a Three-Dimensional  Flow  Problem  ol  Axial-Flow  Turbo- 
machinery," ASME  Transactions,  Vol . 75,  July  19c3,  pp.  789-803. 

8.  Shepherd,  D.  G.,  Principles  of  Turhomachinery , The  MacMillan  Company, 
New  York,  1956. 

9.  Giamati , C.  C.,  Jr.,  Finger,  H.  B.,  "Design  Velocity  Distribution  in 
Meridional  Plane,"  NACA  Special  Publication  36,  196*5,  pp.  255-278. 

10.  Mereiros,  A.  A.,  Hood,  B.  J.,  "Chart  Procedures  for  Design  Velocity 
Distribution,"  NACA  Special  Publication  36,  1965,  pp.  279-295. 

11.  Herzig,  H.  Z.,  Hansen,  A.  G.,  "Three-Dimensional  Compressor  Flow 
Theory  and  Real  Flow  Effects,"  NACA  Special  Publication  36,  1965, 
pp.  365-384. 

Smith,  L.  H.  Jr.,  "The  Radial  Equilibrium  Equation  of  Turbomachinery," 
ASME  Jour na 1 of  Engineering  for  Power , J a nua ry  1966,  pp.  1-12. 

230 


12. 


13. 

14. 

15. 

16. 

17. 

18. 

19. 

20. 

21. 

22. 

23. 

24. 

25. 

26. 


Hawthorne,  W.  R.,  and  Horlock,  J.  H.,  "Actuator  Disk  Theory  of  the 
Incompressible  Flow  in  Axial  Compressors,"  Procedure  of  the  Institu- 
tion of  Mechanical  Engineers,  Volume  176,  Number  30,  1962. 

Hawthorne,  W.  R.,  and  Ringrose,  J.,  "Actuator  Disk  Theory  of  the  Com- 
pressible Flow  in  Free-Vortex  Turbomachinery,"  Thermodynamics  and 
Fluid  Mechanics  Convention  of  Institution  of  Mechanical  Engineers, 
Cambridge  9-10  April  1964. 

Marble,  Frank  E.,  "Three-Dimensional  Flow  in  Turbomachines,"  Section  C 
of  Aerodynamics  of  Turbines  and  Compressors,  Princeton  University 
Press,  1964. 

Dodge,  P.  R.,  "Three-Dimensional  Heat  Transfer  Analysis  Program," 
AiResearch  Manufacturing  Company,  AFAPL-TR-77-64 . 

Binder,  R.  C.,  Advanced  Fluid  Dynamics  and  Fluid  Machinery,  Prentice- 
Hall,  Incorporated , New  York,  1951. 

Velkoff,  Henry  R., "Notes  on  Elements  of  Turbomachinery ," The  Ohio  State 
University,  1969. 

Engelman,  H.  W. , "Unpublished  Lecture  Notes  on  'Advanced  Principles 
of  Energy  Conversion  in  Turbomachinery,  Department  of  Mechanical 
Engineering  of  the  Ohio  State  University,  Spring  Quarter,  1969. 

Bogus,  A.  S.,  "The  Problem  of  Vortex  Flow  in  an  Annulus," 

M,  S.  Thesis,  The  Ohio  State  University,  1970. 

Due,  H.  F.,  Rogo,  C.,  "Advanced  Small  Axial  Turbine  Technology  Program, 
Phase  II  Test  Report,"  Contract  DAAJ02-7 2-C-0117,  February  1974. 

Due,  H.  F.,  Rogo,  C.,  Kosier,  C.  L.,  Jasas,  G.  B.,  "Advanced  Small 
Axial  Turbine  Technology,"  Final  Report  USAAMRDL-TR-77-1 , 

Contract  DAAJ02-72-C-01 1 7 , May  1977. 

Rayle,  R.  E.,  "An  Investigation  of  the  Influence  of  Orifice  Geometry 
on  Static  Pressure  Measurements,"  MIT  S.  M.  Thesis,  Department 
Mechanical  Engineering,  1949. 

Dean,  Robert  C.  Jr.  et  al , Aerodynamic  Measurements,  Gas  Turbine 
Laboratory,  Massachusetts  Institute  of  Technology,  1953. 

Bammert,  K.  and  Sandstede,  H.,  Influences  of  Manufacturing  Tolerances 
and  Surface  Roughness  of  Blades  on  the  Performance  of  Turbines, 

Journal  of  Engineering  for  Power  - Transactions  of  the  ASME, 

January  1976,  pp.  29-36. 


Wysong,  R.,  et  al , "Turbine  Design  System,"  General 
AFAPL-TR-78-92. 


Electric, 


231 


27. 


Keenan,  Joseph  H.,  and  Kaye,  Joseph,  Gas  Tables,  John  Wiley  & Sons, 
Incorporated,  1948. 

28.  John,  James  E.  A.,  Gas  Dynamics,  Allyn  and  Bacon,  1969. 

29.  United  Sensor  and  Control  Corporation,  85  School  Street,  Watertown, 
Massachusetts,  Calibration  Data  - Bulletin  W,  WT 

30.  Keast,  F.  H.  "High  Speed  Cascade  Testing  Techniques,"  Transactions  of 
the  ASME,  Volume  74,  Number  5,  Page  685,  July  1952. 


232 

*U.$. Government  Printing  Office:  1979  — 657-002/16 


